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^A  spectral  infrared  radiative  transfer  model  based  on  the  discrete-ordinate 
method  has  been  developed  for  atmospheres  containing  clouds  and  absorbing  gases 
above,  below  and  within  the  cloud.  The  transfer  model  takes  into  consideration 
the  inhomogeneity  of  the  cloudy  atmosphere,  the  gaseous  absorption  in  scattering 
cloud  layers  and  the  wavenumber  dependence  of  radiative  transfer.  In  addition, 
the  cloud  layer  is  further  divided  into  a number  of  sub-layers  to  account  fdr  the 
non- isothermal  and  inhomogeneous  cloud  characteristics. 

The  spectral  infrared  transfer  program  is  applied  to  VTPR  channels  of  the  — 
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20.  (cont.) 

NOAA  4 satellite  to  simulate  upward  radiances  in  cirrus  cloud  conditions.  Single- 
scattering properties  for  ice  crystals  are  calculated  assuming  ice  cylinders  200 
and  60  unyin  length  and  width,  respectively,  randomly  oriented  in  a horizontal 
plane.  Comparisons  between  satellite  observed  ard  theoretically  simulated  t.rward 
radiances  are  carried  out  for  selected  cirrus  cloud  cases.  Incorporating 
atmospheric  profiles  obtained  from  radiosonde  and  the  observed  cloud  information 
into  the  spectral  transfer  program,  we  show  a systematic  agreement  between  ob- 
served and  computed  upward  radiances.  Systematic  reduction  patterns  of  the  up- 
ward radiance  caused  by  the  increase  of  the  cloud  ice  content  are  demonstrated  for 
VTPR  channels  employing  tropical  and  midlatitude  atmospheric  profiles.  Having 
the  quantitative  relationships  between  upward  radiance  and  ice  contents,  proce- 
dures are  described  for  the  inference  of  the  cloud  ice  content  and  cloud  amount. 
The  proposed  method  has  been  successfully  applied  to  the  three  cirrus  cloud  cases. 

Mdreover,  radiative  transfer  model  of  spectral  infrared  radiation  in  cloudy  A 
atmospheres  is  applied  to  the  most  complete  set  of  radiance  observations  cur- 
rently available  from  the  Nimbus  VI  H I RS  instrument/.  We  show  that  the  effect  of 
the  atmospheric  temperature  profile  is  minimized  by  dividing  the  clear  column 
radiance  to  the  cloudy  radiance.  A method  for  the  estimate  of  cloud  compositions 
from  mul tispectral  H I RS  radiances,  including  both  Short-wave  and  long-wave  COp 
channels,  is  developed.  Theoretical  calculations  bf  the  upwelling  radiance 
at  satellite  altitude  for  a number  of  thicknesses  involving  middle  level  and  cir- 
rus clouds  are  carried  out.  The  resulting  theoretical  radiances  are  then 
parameterized  and  an  empirical  method  to  determine  cloud  type  and  ice  and  water 
content  of  the  cloudy  is  described.  Satellite  passes  from  five  days  over  the 
western  United  States  are  employed  to  test  the  empirical  parameterizations  of  the 
theoretical  results.  Comparisons  of  the  cloud  ty|)e  determinations  with  NOAA  4 
Mosaic  are  shown  to  agree  reasonably  well.  Cloud  ice  and  water  content  derived 
from  the  parameterizations  of  the  HIRS  data  is  compared  with  that  obtained  from 
the  Air  Force  Three  Dimensional  Nephanalysis  (30NEPH)  program.  Examples  of  the 
application  of  this  technique  to  global  mapping  of  ice  and  water  content  are 
displayed. 

In  the  appendix,  a retrieval  technique  is  presented  for  the  determination  of 
the  surface  temperature,  the  thickness  and  transmissivity  of  cirrus  clouds,  and 
the  fraction  of  the  cirrus  cloudiness  by  means  of  four  observed  upwelling  ra- 
diances in  the  10  pm  window  region:  On  the  basis  of  radiative  transfer  calcula- 
tions for  mean  wavenumbers  of  900/  950,  1100  and  1150  cm~'  , assumptions  are  made 
in  the  theoretical  retrieval  analyses  that  water  vapor  effects  above  cirrus 
clouds  are  negligible  and  that  ratios  of  the  transmissivities  are  linear  functions 
of  the  cloud  thickness  Error  analyses  employing  climatological  data  reveal  that 
indepeident  random  errors  in  temperature  and  humidity  profiles  introduce  insigni- 
ficant errors  in  the  four  resulting  parameters.  The  resulting  errors  caused  by 
random  errors  in  the  expected  upwelling  radiances,  however,  depend  upon  their 
standard  deviations.  One?  the  thickness  and  the  transmissivity  at  a given  wave- 
number  of  a irrus  cloud  have  been  determined,  we  illustrate  that  the  vertical 
ice  content  may  be  estimated  assuming  that  ice  particle:  are  randomly  oriented 
in  a horizontal  plane. 
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SECTION  I 


INTRODUCTORY  REMARKS 

The  structure  and  composition  of  clouds  obtained  from  passive 
remote  sensing  from  satellites  on  a routine  basis  have  been  extremely 
limited.  Many  of  the  methods  utilized  are  based  strictly  on  the  statisti- 
cal analysis  of  the  observed  data.  This  may  be  due  to  the  complexity 
of  cloud  interaction  with  the  radiation  field  of  the  atmosphere.  This  is 
especially  evident  for  the  globally  distributed  high,  semi-transparent 
cirrus  clouds.  The  objectives  of  the  satellite  sensing  and  radiative  trans 
fer  project  initiated  in  the  Department  of  Meteorology,  University  of  Utah 
and  sponsored  by  the  Air  Force  Geophysics  Laboratory  are  to  develop  (1)  a 
computer  program  to  simulate  the  transfer  of  infrared  radiation  through 
cirrus  cloudy  atmospheres,  (2)  a passive  satellite  retrieval  technique  to 
estimate  the  vertical  ice  content  profile  of  cirrus  and  (3)  a possible 
three-dimensional  cloud  scene  method  from  a combination  of  available  satel- 
lite sounding  instruments  such  as  those  on  board  Nimbus  VI  and  NOAA  IV. 

In  this  final  report  we  present  some  of  the  highlights  of  results 
obtained  during  the  past  three  years.  Section  2 consists  of  discussions 
on  the  theoretical  foundation  for  the  calculations  of  the  transfer  of 
spectral  infrared  radiation  in  cloudy  atmospheres.  The  infrared  transfer 
model  allows  the  inhomogeneity  of  cloudy  atmospheres  to  be  treated  approxi- 
mately and  it  takes  into  account  the  gaseous  absorption  in  scattering  cloud 
layers  and  the  wavenumber  dependence  of  radiative  transfer.  Application  of 
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the  infrared  transfer  model  to  VTPR  (Vertical  Temperature  Profile  Radio- 
meter) channels  of  the  NOAA  IV  satellite  is  described  in  Section  3. 
Radiative  properties  of  cirrus  clouds  are  investigated  and  possible 
cirrus  cloud  scenes  are  explored.  In  Section  4,  we  apply  the  infrared 
transfer  program  to  HI RS (H i gh  Resolution  Infrared  Sounder)  channels  of 
the  Nimbus  VI  satellite  with  several  combinations  of  middle  and  cirrus 
clouds.  The  theoretical  results  are  parameterized  and  are  then  used  to 
infer  cloud  characteristics  from  actual  H IRS  data.  In  the  appendix, 
we  attach  a paper  (Liou,  1977)  in  conjunction  with  the  development  of  a 
retrieval  technique  for  the  inference  of  cirrus  cloud  characteristics 
based  on  a set  of  synthetic  radiance  calculations  in  the  window  region. 
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SECTION  2 


SPECTRAL  INFRARED  RADIATION  PROGRAM  FOR 
CLOUDY  ATMOSPHERES 

2 . 1 Basic  Equations  and  Approximations 

The  infrared  radiation  program  begins  with  solving  the  transfer 
equation  for  a plane-parallel  cloud  layer  consisting  of  absorbing  gases 
in  local  thermodynamic  equilibrium.  The  basic  equation  describing  the 
monochromatic  infrared  radiation  field  is  given  by 


d I ( T , VI  ) 


= I (t,m)  —7t-  f P (w,u'  )I  (x,u')dvi' 


v — j i vt.ii/ 

at  v c.  _i  v v 


- (l-av)Bv[T(T)], 


(2.1) 


where 


U>  = B /(B  +fi  +nk  ), 
v s.v  s.v  a ,v  v 


B (T)  = 2hc2v3/(ehCv/KT-l), 


(2.2) 


(2.3) 


and  I v represents  the  monochromatic  radiance  of  wavenumber  v.u  the  cosine 

of  the  emergent  angle  with  respect  to  the  zenith,  x the  optical  depth,  Pv 

the  normalized  axially  symmetrical  phase  function,  T the  cloud  temperature 

which  is  a function  of  height  or  optical  depth,  u>v  the  single  scattering 

albedo,  By(T)  the  Planck  function,  h and  K the  Planck's  and  Boltzmann's 

constants,  respectively,  c the  velocity  of  liqht,  and  b the  volume 

s,v  a,v 
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scattering  and  absorption  cross  sections  for  cloud  particles  of  wavenumber 
v,  n the  number  density  of  the  absorbing  gases  within  the  cloud  layer,  and 
kv  the  absorption  coefficient  of  the  gases. 

The  normalized  phase  function  may  be  expanded  into  Legendre 
polynomials  consisting  of  a finite  number  of  terms.  Upon  replacing  the 
integration  in  Eg.  (2.1)  by  summation  according  to  the  Gauss'  quadrature 
formula,  a set  of  first-order  inhomogeneous  differential  equations  are 
derived.  Cy  seeking  the  homogeneous  and  particular  solutions  of  the 
differential  equations  as  described  by  Chandrasekhar  (1950),  the  complete 
solutions  of  the  scattered  radiance  for  a given  discrete-stream  j assuming 
an  isothermal  cloud  temperature  Tc  may  be  written  (Liou,  1 973 ; 1 974 ) 

-k  t 


li.)  = I L $ ( u • ) e 
i m mm  l 


* W- 


(2.4) 


where  r denotes  summation  over  the  2n  discrete  streams  employed,  <j>m  and 
km  are  the  eigenfunction  and  eigenvalue  of  the  differential  equations 
whose  values  depend  upon  the  phase  function  and  single-scattering  albedo, 
and  Ln)  are  a set  of  constants  of  proportionality  to  be  determined  from 
the  radiation  boundary  conditions  above  and  below  the  cloud  layer. 

The  upward  and  downward  radiances  arising  from  the  molecular 
absorption  and  emission  reaching  the  cloud  bottom  and  top,  respectively, 
can  be  obtained  by  solving  the  transfer  equation  for  a non-scattering 
atmosphere  in  local  thermodynamic  equilibrium.  They  are  given  by 

= 1 Bv[T(z)]d  Tv(z,zt;-ui ),  (2.5) 
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I V ,u.)  = B (T  )T  (2  ,0;p.) 

v D 1 v S v D 1 
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+ / B [T(z)]d  T (z, ,z;u.),  (2.6) 
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where  T$  is  the  surface  temperature,  and  zt  and  z^  are  cloud  top  and  base 
heights,  respectively,  and  the  monochromatic  transmittance  is 

1 22 

T (z-.z.iy.)  = exp  [-— / k (z)n(z)dz].  (2.7) 

v 2 1 i y .•  , v 

1 Z1 


Assume  that  the  variation  of  the  Planck  function  with  respect  to  the 
wavenumber  is  much  smaller  than  that  of  the  transmission  function  T . 
Thus,  upon  multiplying  Eqs.  (2.5)  and  (2.6)  by  the  instrumental  slit 
function  $(v)  and  performing  the  wavenumber  integration  over  the  spectral 
interval  (v-]^),  we  have 


dv 

Av 


■ < VT(2)]d  TJ2-2t  ••V- 


ow  - 1 k'wm  & 

vi 


(2.8) 


Zb 

+ / B [T(z)]d  T (z.,z;y.),  (2.9) 

0 Av  Av  b 1 

where  the  spectral  transmittance  is  defined  by 
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v2  z2 

T A„(z-.2, in. ) = / exp  [ — — f k (z)n(z)dz]  *(v)  (2.10) 

AV  2 1 1 U,  2 v Av 

We  first  examine  the  transmission  function  (or  transmittance). 

For  a given  height  z,  the  transmittances  are  normally  available  in  the 
form 


v2 


TAv(z)  = TAv(z’“;1)  = 1 exP  [-  ; k (z)n(z)dz]  *(v)  (2.11) 

°° 


In  reference  to  Eqs.  (2.1)  and  (2.2),  we  notice  that  the  gaseous 
absorption  coefficient  is  needed  to  carry  out  the  transfer  of  infrared 
radiation  in  cloud  layers  composed  of  absorbing  gases.  However,  kv  is  known 
only  through  the  absorption  line  parameters  of  gases  and  its  values  vary 
greatly  with  wavenumber  in  a small  spectral  interval.  It  is  very  diffi- 
cult, if  not  impossible,  to  carry  out  line-by-line  calculations  including 
scattering  contributions  of  cloud  particles.  Thus,  a simpler  approach 
for  the  gaseous  absorption  in  a scattering  layer  would  be  to  make  use  of 
the  known  transmittances  which  have  been  obtained  to  a good  accuracy  by 
means  of  line-by-line  calculations  for  inhomogeneous  atmospheres  in  con- 
junction with  satellite  sensing.  Recognizing  the  definition  of  the 
vertical  transmittance  in  Eq.  (2.11),  we  may  approximate  it  by 
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where  the  vertical  path  length 
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k.  may  be  thought  of  as  an  equivalent  absorption  coefficient,  w.  is  the 

J J 

weight  and  M denotes  the  total  number  of  finite  terms  in  the  fitting  of 
the  transmittances.  Once  k-  and  w.  have  been  determined,  we  may  consider 

J J 

the  transfer  of  spectral  infrared  radiation  as  monochromatic  in  the  sub- 
spectral  interval  j and  carry  out  transfer  calculations  in  a cloud  layer 
M times  with  a new  single  scattering  albedo  <3vJ.  Here,  we  assume  that 
the  phase  function  and  the  absorption  and  scattering  cross  sections  of 
ice  crystals  are  independent  of  the  wavenumber  within  a spectral  interval. 
These  assumptions  are  justified  in  view  of  the  relatively  slow  varying 
refractive  indices  of  ice  in  the  infrared  regions  as  evident  in  the  next 
section. 

Without  postulating  any  assumptions  and  approximations,  the  upward 
radiance  at  the  satellite  point  of  view  in  completely  cloudy  conditions  may 
be  expressed  by 
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r i 

where  Ty  denotes  the  monochromatic  transmittance  associated  with  the  up- 
ward radiance  I*  at  the  cloud  top.  The  second  term  on  the  right-hand 
side  is  for  the  clear  atmosphere  above  the  cloud  layer  and  it  is  an  exact 
expression.  However,  the  first  term  requires  approximations  in  order  to 

7 

LL 


i 


I 


i 

i 

i 


take  the  scattering  of  cloud  particles  into  account. 

On  the  basis  of  the  exponential  fit  to  the  transmission  function 
described  in  Eq.  (2.12),  let  ij  (j=l ,2,. . . ,M)  be  the  resulting  radiances 
calculated  from  the  transfer  program  involving  cloud  layers.  Eq.  (2.13) 
can  be  approximated  by 


+ / BA  LT(z)]  dTA  (-.z^),  (2.14) 

2t 


where  at  the  cloud  top  is  to  be  evaluated  from  the  upward  and  downward 
radiances  arising  from  the  molecular  absorption  and  emission  reaching  the 
cloud  bottom  and  top  using  the  equivalent  absorption  coefficients  and 
weights  derived  from  the  exponential  fit  noted  in  Eq.  (2.12),  and  T-  is 

J 

associated  with  the  exponential  function  given  by  Eq.  (2.12).  The  success 
of  carrying  out  such  a semi -monochromatic  infrared  transfer  calculation 
depends  solely  on  the  reliability  of  fitting  the  satellite  transmission 
functions  for  the  entire  atmosphere. 

Since  the  variation  of  the  radiative  property  of  clouds  over  a 
small  spectral  interval  is  much  smaller  than  that  of  the  transmittance  and 
since  the  spectral  transmi ttances  for  satellite  channels, TAv,  are  available 
from  the  top  of  the  cloud  to  the  top  of  the  atmosphere,  we  postulate  from 


Eq.  (2.13)  that 


* ‘L'vi’  V"V“(> 


(2.15) 
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where 


(2.16) 


and  semi -monochromatic  upward  radiances  it  are  to  be  evaluated  from  the 

J 

procedures  described  below.  The  approach  described  here  makes  use  of 
the  known  spectral  transmission  functions  above  and  below  the  cloud  and 
utilizes  the  exponential  fitting  program  for  the  cloud  layer.  If  the 
cloud  layer  is  a black  body,  then  Eq.  (2.15)  is  an  exact  expression  where 
1^  will  be  given  by  the  Planck  function  of  the  cloud  top  temperature. 

For  moderately  thick  and  thick  clouds,  we  would  think  that  the  approxi- 
mation is  a good  one  since  cloud  particles,  not  the  gases  within  the 
cloud,  dominate  the  transfer  processes. 

The  upward  radiance  at  the  satellite  point  of  view  in  clear 
conditions  is  simply 


+ / BJ TU>] 


(2.17) 


Thus,  if  within  the  f ield-of-view  of  the  satellite  radiometer  there  is  n 
portion  of  cloudiness,  the  upward  radiance  at  the  top  of  a partly  cloudy 


atmosphere  is  then  given  by 


*Av  ("*»,i)  = nIAv  + (1*n)IA\,NC(~’Wi)- 


(2.18) 


2.2  Radiative  Transfer  in  Non- Isothermal  and  Inhomogeneous  Cloud 

Layers 

The  solution  of  the  infrared  radiative  transfer  equation  given 
by  Eq.  (2.4)  is  applicable  only  to  isothermal  and  homogeneous  cloud  layers. 
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The  question  of  inhomogeneity  of  cirrus  cloud  compositions  is  a difficult 
one.  In  the  first  place,  there  have  not  been  many  observations  available. 
For  the  purpose  of  radiative  transfer  calculations,  it  seems  that  averaged 
properties  of  cloud  compositions  may  be  appropriate.  However,  the 
effect  of  the  non- isothermal  structure  of  clouds  on  their  radiative  pro- 
perties is  more  critical.  Physically,  a cold  cloud  top  would  reduce  the 
transmitted  radiances.  It  is  conceivable,  therefore,  that  an  overest inat ion  of 
upward  radiances  is  likely  to  take  place.  Although  our  main  concern  is 
the  non- isothermal  structure,  the  transfer  program  described  below  may 
be  utilized  to  investigate,  if  desirable,  the  inhomogeneous  properties 
of  clouds  as  wel 1 . 

In  •'eference  to  Figure  2.1,  the  cirrus  cloud  layer  is  divided 

into  a number  of  sub-layers  each  of  which  is  considered  to  be  isothermal 

and  homogeneous . The  optical  depth  is  evaluated  from  the  cloud  top  to 

the  bottom  of  the  sub-layer.  The  index  « is  used  to  denote  the  number  of 

sub-layer.  We  now  apply  the  solution  of  the  radiative  transfer  equation 

given  by  Eq.  (2.4)  to  each  sub-layer  and  sub-spectral  interval  j to 

obtain  (>:  denotes  summation  over  discrete-streams,  (-n,n)) 
m 


m (u.)e 
m mm  i 


. a 

-k  r 
m 


+ B, 


<r>. 


In  order  to  determine  the  unknown  coefficients  L^,  the  radiation 
continuity  relationships  are  needed  in  addition  to  two  radiation  boundary 
conditions  specified  in  Eqs.  (2.8)  and  (2.9).  At  the  cloud  top  the  down- 
ward radiance  has  to  be  equal  to  that  from  the  molecular  atmosphere  above, 
so  that 
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Figure  2.1.  Radiative  transfer  model  for  an  inhomogeneous  cirrus  cloudy 
atmosphere.  The  cirrus  cloud  layer  is  divided  into  several 
sub-layers  each  of  which  is  considered  to  be  isothermal  and 
homogeneous. 


L 


n 


i 


h 1 > = rAv+(z  *_lJi)wj  • 


(2.20) 


Between  the  layers,  the  radiances  from  all  directions  must  be  continuous. 


Thus, 


Ij  = Ij  f 1 ),a=l  ,2,. . . ;N-1 , 


where  N denotes  the  total  number  of  sub-layers.  Lastly,  the  upward 
radiance  has  to  be  equal  to  that  from  the  molecular  atmosphere  below 


(2.21) 


to  give 


Ij  - Iav,(^.*ui)"j  • 


(2.22) 


Upon  inserting  the  radiance  solution  expressed  by  Eq.  (2.19)  into 


Eqs.  (2.20)  - (2.22),  we  have 


■ 1 /(W  - Bj  <Tc>- 


t+1  „£+l 


s UKM  + L;  45  (wi)]  ■ Bj  (Tj  ')  - Bj  (Tj), 


£LM(^i)-1  /(W-Bj  <Tc>’ 


(2.23) 


(2.24) 


(2.25) 


where 


-kV 


£+1  £ 

.£+1,  . £,  . "km  T 

fim  = ‘ *n>i>e 


(2.26) 


(2.27) 


0' 


Eqs.  (2.23)  - (2.25)  represent  a system  of  2n  x N linear  equations  from 
which  the  L((|  coefficients  may  be  determined  using  a matrix  inversion 
technique.  We  may  rewrite  Eqs.  (2.23)  - (2.25)  in  a compact  matrix  form 

* L * B,  (2.28) 


where  the  unknown  coefficients  of  proportionality 


(2.29) 


The  matrix  denoting  the  contribution  due  to  cloud  emission  and  upward 
and  downward  radiances  reaching  the  cloud  base  and  top,  respectively. 


1 . 
J 

*“fV 

-BJ 

(Tj) 

■j 

- 8. 

J 

(tJ) 

BJ 

(7c>  - Bj 

<tJ> 

1 

1 

BJ 

(Tc>  - Bj 

(tJ) 

| 

'j 

^ V+l‘l  * 

•BJ 

a?)  i 

c i 

’j 

- B. 

J 

(T^)  i 

(2.30) 
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and  the  2n  x N by  2n  x N matrix 


[ 


*n(-l,n} 

’ ^{"ul) 

>-i('“n) 

••  Yi('v‘n)5-n('v‘n) 

..  .*( 

Y ](p  ) 

-n  n 

. . Y^ ( P )5  2(p  ) 
n n -n  n 

..  «2( 
n 

M ) 

n 

*-nS> 

••  *!<*!> 

1 

••  %<“„> 

4 


(2.31) 


In  Eq.  (2.31),  the  blank  spaces  denote  zero  elements.  A similar  procedure 
has  been  employed  by  Liou  (1975)  for  the  transfer  of  solar  radiation  in 
inhomogeneous  atmospheres.  Once  have  been  determined,  they  can  be 
inserted  in  Eq.  (2.19)  to  obtain  the  sub-spectral  radiance  distribution 
within  each  sub-layer. 


2.3  Discussions  on  Numerical  Procedures 

Calculations  for  single  scattering  parameters  were  first  carried 
out  for  ice  cylinders  randomly  oriented  in  a horizontal  plane  according  to 
a previous  theoretical  model  for  ice  crystal  clouds  (Liou,  1972).  Based 
on  observations  of  Weickmann  (1949)  and  experimental  results  from  Heymsfield 
and  Knollenberg  (1972),  a mean  length  of  200  pm,  a mean  radius  of  30  pm  and 
a mean  concentration  of  0.05  cm'^  along  with  the  corresponding  refractive  in- 
dices for  ice  (Schaaf  and  Williams,  1973)  were  assumed  in  the  scattering  cal- 
culations. The  single  scattering  parameters  for  the  water  cloud  utilized  a 
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drop  size  distribution  developed  by  Feddes  and  Smith  (1974).  This  distri- 
bution is  an  exponential  fitting  of  the  best  published  measurements  taken 
from  cloud  physics  literature.  Although  this  fittinq  was  done  for  ten  cloud 
types,  only  the  altostratus  cloud  distribution  is  used.  In  the  calcula- 
tion of  the  single  scattering  parameters  for  altostratus  clouds,  a liquid 

_ O 

water  content  of  0.15  gm  m was  used.  The  real  and  imaginary  parts  of 
the  refractive  index  for  water  were  taken  from  Hale  and  Querry  (1973). 
Single-scattering  computations  were  carried  out  for  the  central  wavenumber 
of  each  channel . 

The  phase  functions  are  expanded  in  a series  of  Legendre  poly- 
nomials in  the  form 
N 

P (u.u')  ~ £ tiPJu)P.(v'),  (2.32) 

1=0  *’  x * 

where  are  obtained  from  the  orthogonal  properties  of  the  Legendre 
polynomials.  The  expanded  form  is  to  be  used  in  the  transfer  calcula- 
tions. 

To  derive  a set  of  equivalent  absorption  coefficients  for 
use  in  cloud  scattering  calculations,  an  exponential  fit  of  the  form 
denoted  in  Eq.  (2.12)  was  applied  to  transmission  curves  of  the  CO^ 
and  HpO  channels.  Note  that  the  transmission  curves  take  into  account 
absorption  by  H20,  C02  and  03  simultaneously,  with  each  gas  distributed 
differently,  and  the  instrumental  slit  function.  These  curves  cover 
the  entire  atmosphere  from  0.01  mb  to  the  surface,  which  introduces 
a significant  pressure  variation. 

A simplified  fitting  routine  was  developed  which  guaranteed  the 
required  restriction  on  Wj  and  kj  (both  have  to  be  positive).  The  numerical 
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scheme  is  based  on  the  followinq  iterative  process.  An  initial  T (u.) 

Av  I 

and  Wj  are  chosen  (initially  Wj  = 1 ) and  kj  is  generated  from  (see  Eq. 

(2.12)) 


«.n 


(2.33) 


-k^u 

The  curve  T^v(u)  = e is  then  generated  point  by  point;  as  each 
new  point  is  generated  it  is  compared  with  the  true  value  TAv(u),  which  is 
the  known  transmi ttance.  to  make  sure  it  meets  the  following  conditions: 


T!  (u)  < T (u) 

Av  Av 


dT*  (u)  dT  (u) 

Av av 


du 


du 


J 


(2.34) 


The  second  condition  requires  only  a rough  approximation  for  the  derivative. 
If  either  of  these  are  not  met,  a new  w^  is  generated,  (i.e.,  w1  = w^  x 
constant)  so  that  the  new  w1  is  slightly  less  than  the  old  one,  and  the 
above  is  repeated.  If  no  w-j  can  be  found,  a new  Ui  is  selected  (larger 
than  the  old  point).  When  a good  w^  and  k-j  are  found,  the  above  process 
is  then  repeated  to  produce  a w^  and  k?  using  values  of  TAv(u)  - TAv(u) 
instead  of  TAv(u).  The  number  of  points  needed  to  fit  the  curve  then 
depends  only  on  the  size  of  the  error  tolerance  used  and  on  the  range  of 
u over  which  TAv(u)  is  to  be  fitted.  The  exponential  fitting  approximation 
is  restricted  in  the  scattering  cloud  layer,  so  that  the  gaseous  absorption 
within  the  cloud  layer  may  be  included.  Note  here  that  spectral  transmit- 
tances  are  available  both  above  and  below  the  cloud. 

Once  the  scattering  and  absorption  properties  of  ice  crystals 
and  the  spectral  absorption  characteristics  of  gases  have  been  determined. 
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we  may  evaluate  the  volume  single-scattering  properties  for  a mixture  of  ice 
crystals  and  gases  within  the  clouds.  According  to  Eq.  (2.2)  and  the  expo- 
nential fitting  program  described  above,  the  s ingle-scattering  albedo,  in- 
cluding gaseous  contribution  for  a spectral  interval  whose  central  wavenumber 
is  v,  may  be  written 

*1  - Ss  /(be+nkj) , j = 1,  . . . ,N . (2.35) 

where  t>e  represents  the  volume  extinction  cross  section.  The  total  optical 
depth  for  a cloud  with  a thickness  of  az  is  given  by 

= 6eAz  + ukj,  j = 1 N,  (2.36) 

where  we  omit  the  subscript  v in  the  right-hand  side  of  these  two  equations. 
Note  that  u is  the  gaseous  pathlenqth.  The  amount  of  i^O  within  clouds  is 
estimated  by  using  a mean  temperature  and  assuming  saturated  conditions.  CO., 
concentration  is  assumed  uniformly  distributed  in  the  atmosphere.  Thus,  the 
pathlength  of  gases  within  the  cloud  may  be  obtained.  Finally,  it  is  also 
assumed  that  the  spectral  phase  function  is 

p(j(u,w' ) = pv(n.u ' ) * (2.37) 

as  described  previously.  The  single-scattering  albedo,  the  total  optical 
depth  for  a given  thickness  az  and  the  phase  function  in  the  form  of  Legendre 
polynomials  are  the  basic  parameters  for  carrying  out  the  spectral  radiative 
transfer  calculations. 

The  monochromatic  transfer  program  first  generates  the  required 
eigenvalues  for  a set  of  homogeneous  differential  equations  from  the 
single-scattering  albedo  and  the  expanded  phase  function.  Liou  (1973) 
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pointed  out  the  mathematical  and  numerical  ambiguities  involving 
Chandrasekhar ' s (1950)  method  for  searching  the  eigenvalues  correspond- 
ing to  the  2n  associated  homogeneous  differential  equations,  and  developed 
a matrix  method  directly  from  the  associated  homogeneous  differential 
equations  for  the  eigenvalue  problem.  Recognizing  the  symmetric  relation- 
ships between  the  matrix  elements,  Asano  (1975)  further  developed  the  ana- 
lytic procedures  for  the  reduction  of  the  rank  of  the  matrix.  These 
developments  thus  allow  the  accuracies  of  eigenvalues  to  be  improved  and 
the  computer  time  involved  greatly  reduced.  Computations  presented  in 
this  paper  follow  the  numerical  procedures  outlined  by  Liou  (1973)  and 
Asano  (1975)  for  eigenvalue  problems  in  the  discrete-ordinate  method  for 
radiative  transfer. 

Secondly,  the  Drogram  calculates  the  clear  column  angular  upward 
and  downward  radiances  arising  from  the  absorption  and  emission  of  gases 
(in  forms  of  transmi ttances)  above  and  below  the  cloud  as  well  as  the 
surface  contribution.  Each  model  atmospheric  profile  is  divided  into  101 
levels  from  the  surface  to  100  km  by  means  of  linear  interpolation  between 
the  available  data  on  pressure,  temperature  and  humidity.  The  non-isother- 
mal  cloud  temperatures  are  assumed  the  same  as  those  of  the  surroundings. 

The  program  then  solves  the  unknown  coefficients  of  proportion- 
ality of  radiance  solutions  for  an  inhomogeneous  system  of  equations 
described  in  Section  2.2.  After  proper  summation  over  the  sub-spectral 
weight  derived  from  the  exponential  fit,  the  transmitted  and  reflected 
spectral  radiances  at  the  cloud  top  and  bottom,  respectively,  can  be 
obtained.  Since  the  infrared  computer  program  emoloys  16  discrete  streams 
in  the  calculations,  interpolations  of  the  16  radiance  values  are  finally 
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needed  to  obtain  0°,  180°  and  1 i mb  (90°)  directions.  Note  that  0°  and 
180°  directions  are  used  to  define  the  transmissivity  and  reflectivity 
of  clouds. 


r 


SECTION  3 


APPLICATION  TO  NOAA  IV  VTPR  CHANNELS 

3. 1 Characteristics  of  VTPR  Channels 

In  this  section  we  investigate  effects  of  cirrus  clouds  on 
NOAA  IV  VTPR  (Vertical  Temperature  Profile  Radiometer)  channels  and  explore 
cirrus  cloud  scenes  utilizing  upwelling  radiance  observations  from  these 
channels. 

The  NOAA  4 VTPR  instruments  consist  of  six  channels  in  the  15  ym 
band  of  C02>  one  channel  in  the  window  region  at  12  ym  and  one  channel  in 
the  rotational  band  of  water  vapor.  Table  3.1  shows  the  nominal  character- 
istics of  the  eiqht  filters  in  the  VTPR  instruments.  The  center  wavenumbers 
arp  employed  to  calculate  the  single  scattering  properties 


Table  3.1.  Nominal  spectral  intervals  for  VTPR  channels, 
after  McMillin  et  al . (1973). 


Channel 

Center  wavelength 

Half-width 

(uni) 

(cm  ]) 

(cm-1 ) 

1 

14.96 

668.5 

3.5 

2 

14.77 

677.5 

10 

3 

14.38 

695.0 

10 

4 

14.12 

708.0 

10 

5 

13.79 

725.0 

10 

6 

13.38 

747.0 

10 

7 

18.69 

535.0 

18 

8 

11.97 

833.0 

10 
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of  ice  crystals  which  are  presumably  valid  for  the  entire  spectral  region. 
Atmospheric  transmi ttances  for  the  eight  spectral  intervals  are  available 
for  the  zenith  angles  of  0°  and  23°47'  (McMillin  et  al.,  1973).  Accord- 
ing to  their  report,  transmi ttances  for  C02  are  based  upon  calculations 
made  by  Orayson  (1971)  employing  the  point-by-point  method  for  a number 
of  temperature  profiles.  From  these  transmi ttances,  transmi ttances  for  a 
given  temperature  profile  may  be  obtained  by  means  of  interpolation.  Ozone 
transmittances,  which  are  a minor  correction,  are  calculated  from  a line-by- 
line technique  for  the  lower  atmosphere  using  climatological  ozone  profiles 
(McClatchey  et  al.,  1971).  As  for  water  vapor,  a procedure  proposed  by  Weinreb 
and  Neuendorffer  (1973)  is  utilized  to  obtain  the  transmittances  of  selec- 
tive absorption  in  inhomogeneous  atmospheres.  The  effect  of  the  water 
vapor  continuum  is  also  included  based  on  laboratory  measurements  (Bignell, 
1970).  Thus,  the  transmittance  of  a VTPR  channel  for  a given  height  z is 
given  by 


C0„  H90  0, 

T (z)  = T 2(z)T  2 (z)T  3(z) 
Av  Av  Av  Av 


(3.1) 


Figure  3.1  depicts  the  weighting  functions  of  the  VTPR  channels. 
Since  clouds  are  normally  tropospheric  in  origin  and  are  located  below  about 
200  mb  or  so,  we  would  anticipate  that  their  effects  on  upward  radiances  ob- 
served in  channels  1 , 2 and  3 of  the  C02  band  are  relatively  insignificant. 

There  are  only  two  sets  of  transmittances  being  generated  for 
VTPR  instruments  corresponding  to  0°  (nadir)  and  23°47'  zenith  angles. 
However,  spectral  radiative  transfer  for  cirrus  clouds  require  the  clear 
column  angular  upward  and  downward  radiances  reaching  the  cloud  base  and 
top,  respectively.  To  evaluate  these  radiances,  the  angular  dependent 
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transmittances  from  0°  to  90°  zenith  angles  are  needed.  The  maximum 
VTPR  scan  angle  is  about  40°  zenith  angle.  Thus,  angular  transmittances 
other  than  0°  and  23°47‘  are  also  needed  for  the  interpretation  of  the  ob- 
served upward  radiances.  According  to  McMillian  (1976,  personal  communi- 
cation), the  angular  transmittances  for  the  15  pm  CC>2  channels  are  approx- 
imately given  by 


tAv(z»w)  = [TAxl(z)] 


1 


Av' 


(3.2) 


where 

1 0J  f2 

f^  = [sec(cos  p)]  . (3.3) 

The  empirical  coefficient  f2  was  derived  from  a four  term  polynomial 
fitting.  It  is  related  to  the  pressure  and  temperature  for  a given  level 
and  certain  empirically  determined  constants.  For  water  vapor  channels  7 
and  8,  however,  angular  dependent  transmittances  are  approximated  by 

T (z,u)  = [Ta  (z)]sec(0‘77  cos  \) . (3.4) 

Av  Av 

Verification  of  the  empirical  angular  transmittances  was  done  by 
comparing  the  calculated  and  observed  radiances  for  various  scan  angles. 

It  appears  that  up  to  about  40°,  they  are  reasonably  reliable.  Since  no 
additional  information  is  available,  we  have  used  these  empirical  equations 
to  evaluate  the  angular  transmittances  from  0°  to  90°  zenith  angles  for  the 
eight  VTPR  channels.  It  is  anticipated  that  errors  may  be  produced  for 
large  zenith  angles  close  to  the  limb.  But  we  also  note  that  contributions 
of  radiances  close  to  the  limb  are  insignificant  owing  to  the  very  long 
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gaseous  path  length. 


Table  3.2  lists  optical  properties  for  randomly  oriented  ice  cy- 
linders based  on  a scattering  program  developed  by  Liou  (1972).  The  real  and 
imaginary  parts  of  the  refractive  indices  are  taken  from  the  recent  measure- 
ments by  Schaaf  and  Williams  (1973).  We  see  that  imaginary  parts  decrease 
from  the  shorter  wavelength  to  the  longer  wavelength  in  this  part  of  the  in- 
frared spectrum  with  the  window  wavelength  having  the  largest  value.  However, 
the  window  wavelength  is  also  associated  with  the  smallest  value  of  the  real 
part  of  the  refractive  index.  The  ice  cylinders  appear  to  have  the  largest 
extinction  cross-section  and  largest  single  scattering  albedo  in  the  rotational 
band  wavelength  owing  to  the  smallest  imaginary  part.  The  fact  that  the  ex- 
tinction and  single-scattering  albedo  vary  insignificantly  in  these  wavenumbers 
is  probably  because  of  the  large  particle  size  considered  in  this  study. 

Table  3.2  Optical  properties  of  ice  cylinders  (length  200  vim, 
width  60  urn)  for  the  central  wavenumbers  of  VTPR 
channels. 


Channel 

v(cm_1 ) 

nr 

ni 

r°o 

aext(10~4cn’2) 

1 

668.5 

1.573 

0.178 

0.527 

2.94 

2 

677.5 

1.581 

0.194 

0.528 

2.94 

3 

695.0 

1.583 

0.232 

0.530 

2.94 

4 

708.0 

1.577 

0.262 

0.532 

2.92 

5 

725.0 

1.556 

0.305 

0.533 

2.92 

6 

747.0 

1.522 

0.349 

0.534 

2.92 

7 

535.0 

1.503 

0.074 

0.554 

3.12 

8 

833.0 

1.259 

0.409 

0.531 

2.86 
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3 • 2 Comparison  Between  Satellite  Observations  and  Theoretical 

Simul ations  for  Selected  Cirrus  Cloud  Cases 

Three  cirrus  cloud  cases  under  the  NOAA  satellite  pass  were 
selected  for  comparison  between  the  present  theoretical  calculations 
and  satellite  observations . Figures  3.3a,  3.3b,  and  3.3c  represent 
the  case  studies  to  which  the  cloud  model  was  applied.  Each  case  study 
is  described  to  include  the  atmospheric  profile,  the  synoptic  situation, 
and  the  cloud  composition  determination.  These  cases  were  selected 
from  a group  of  cases  developed  by  the  AFGL  during  1974  and  1975.  These 
three  case  studies  were  used  to  test  the  cloud  model  program  under  three 
conditions  and  include  cirrus  with  low  clouds  (Jan.  11,  1974),  thin 
cirrus  with  no  low  clouds  (Feb.  1,  1974),  and  thick  cirrus  with  no  low 
clouds  (Jan.  22,  1975). 

Each  figure  depicts  the  profile  of  temperature,  pressure,  and 
water  vapor  versus  height.  Each  parameter  was  interpolated  to  101  levels 
used  in  the  clear  column  radiance  calculations.  Input  into  each  profile 
was  the  actual  temperature  and  water  vapor  profiles  of  the  nearest  radio- 
sonde station  to  the  point  of  interest  in  time  and  space  linearly  inter- 
polated to  the  required  pressure  heights.  Above  the  level  of  the  radiosonde 
the  temperature,  water  vapor,  height  and  pressure  profiles  compiled  by 
McClatchey  et  al . (1971)  for  a midlatitude  winter  atmosphere  were  used  for 
each  case  study. 

The  position  and  ice  content  of  the  cloud  for  each  case  were 
estimated  from  airborne  instrumentation  equipped  with  a Knollenberg 
probe  (Knollenberg,  1970)  used  by  AFGL  during  the  cloud  physics  experi- 
ment. These  estimates  were  obtained  by  an  aircraft  spiralling  downward 
from  10  km  at  the  approximate  time  of  the  satellite  pass.  Objective 
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As  in  Figure  3.3(a)  except  for  a cirrus  cloud  case  consisting  of  a high  cirrus 
with  a stratocumulus  below. 


RICHMOND.  INDIANA 
(F«b  I.  1974) 


As  in  Figure  3.3(a)  except  for  a cirrus  cloud  case  consisting  of  multi-layered 
cirrus  clouds. 


analyses  of  ice  contents  were  also  made.  The  ice  content  analyses  involve 
best  estimates  derived  from  detailed  aircraft  observations  consisting  of 
visual  observations,  snow  intensity,  particle  size  replica,  visibility, 
temperature  and  photography  (Conover,  1976,  personal  communication). 

From  all  these  pieces  of  information  as  functions  of  the  altitude,  Conover 
(with  the  consultation  of  AFGL  cloud  physists  headed  by  R.  Cunningham) 
derived  the  total  vertical  ice  contents  used  in  this  investigation. 

To  use  the  observed  cloud  information  in  the  theoretical 
transfer  calculations,  the  vertical  ice  content  for  the  ice  cloud  model 
employed  in  this  study  has  to  be  defined.  For  randomly  oriented  ice  cyl- 
inders in  a horizontal  plane,  it  is  given  by 

IC  = itr^i  pj  N az,  (3.5) 

where  r and  e,  are  the  radius  and  the  length  of  the  cylinder,  respectively, 
Pj  the  density  of  ice,  N the  number  density  and  Az  the  thickness  of  the 
cloud.  As  described  earlier,  the  vertical  ice  content  and  the  cloud 
thickness  are  available  from  aircraft  observations.  Since  the  radius 
and  the  length  of  the  ice  cylinder  have  been  specified  to  obtain  the 
single-scattering  parameters,  it  follows  that  the  "equivalent"  particle 
number  density  can  be  derived.  Once  the  number  density  has  been  obtained, 
the  optical  depth  of  ice  within  the  cirrus  cloud  under  consideration  is 
simply  equal  to  oext  N az.  Together  with  the  optical  depth  of  the  gases 
within  the  cloud,  the  total  optical  depth  represents  one  of  the  fundament- 
al parameters  in  transfer  calculations.  Having  the  atmospheric  temperature 
and  water  vapor  profiles  from  radiosonde,  and  the  cloud  structure  and 
composition  specified,  transfer  calculations  may  now  be  carried  out  for 
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these  cases. 


The  satellite  data  depicted  on  each  figure  is  the  scanning 
radiometer  (SR)  data  representing  the  area  sensed  at  one  scan  position 
of  the  VTPR  sensor,  where  the  zenith  angle  is  that  angle  that  the  center 
of  the  VTPR  scan  position  makes  with  the  earth's  local  zenith.  The  upper 
presentation  is  the  temperature  reading  obtained  from  the  scanning  radio- 
meter, while  the  lower  presentation  is  the  visible  channel  of  the  scanning 
radiometer.  Brief  descriptions  on  each  case  are  in  order. 

January  11,  1974 

The  area  of  the  United  States  east  of  the  Rocky  Mountains  was 
under  the  influence  of  a weak  low,  which  at  500  mb  had  its  center  in 
north  central  Canada.  The  trough  axis  was  very  flat  extending  through 
North  Dakota  southwestward  into  western  Texas.  A weak  surface  low  was 
centered  in  western  Pennsylvania  with  a weak  cold  front  extending  south- 
westward  to  New  Orleans.  A weak  warm  front  extended  due  eastward.  Pre- 
cipitation was  spread  along  the  frontal  system  from  Alabama  to  the  east 
coast.  Chesterfield,  South  Carolina  34.65  N - 80.20  W was  located  under 
the  cirrus  shield  of  the  frontal  system.  The  location  had  an  observation 
of  overcast  stratocumulus  with  the  top  at  2 km  and  a 6/10  cirrus  layer 
from  8.8  to  9.8  km.  The  temperature,  pressure  and  water  vapor  profiles 
were  obtained  from  the  1200Z  Athens,  Georgia  sounding.  The  time  of  the 
observation  and  the  satellite  pass  time  was  1521Z.  The  cirrus  cloud  ice 
content  was  estimated  to  be  11  gm  m"2  in  the  VTPR  field  of  view.  The  SR 
field  of  view  of  the  VTPR  spot  had  a mean  brightness  of  181  and  a mean 
temperature  of  281. 3°K.  These  readings  suggest  that  the  cirrus  is  semi- 
transparent in  both  channels.  Most  of  the  observed  brightness  was  probably 
contributed  by  the  stratocumulus  since  the  cirrus  was  semi-transparent  to 
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the  aircraft  observer. 

F ebruary  1,  1974 

The  area  of  the  United  States  east  of  the  Rocky  Mountains  was 
under  the  influence  of  zonal  flow  at  500  mb  with  a surface  high  pressure 
centered  near  the  Minnesota-Canadian  border.  A weak  frontal  system  was 
oriented  east-west  across  central  Alabama  to  a low  pressure  in  western 
Oklahoma.  A major  area  of  precipitation  was  apparent  over  Michigan  west- 
ward to  North  Dakota.  Richmond,  Indiana  39.75  N - 84.83  W was  located 
under  cirrus  clouds  south  of  the  main  precipitation  area.  The  location 
had  an  observation  of  a thin  7/10  cirrus  cover  from  9-9.8  km.  The  temper- 
ature, pressure  and  water  vapor  profiles  were  obtained  from  the  12007 
Dayton,  Ohio  sounding.  The  time  of  the  observation  and  satellite  pass 
time  was  1526Z.  The  cirrus  cloud  ice  content  was  estimated  to  be  2 gm  nf*- 
in  the  VTPR  field  of  view.  In  this  case,  the  SR  means  are  58  in  the  visible 
channel  and  263. 9°K  in  the  IR  channel.  The  cirrus  is  again  transparent  in 
both  channels  since  the  mean  readings  are  very  close  to  the  surface  con- 
ditions. Note  that  during  the  above  two  days,  the  cloud  physics  obser- 
vations were  under  the  N0AA  2 satellite.  Inspection  of  the  VTPR  trans- 
mittances  of  N0AA  2 and  N0AA  4 indicates  insignificant  differences.  Thus, 
we  have  consistently  used  the  N0AA  4 VTPR  transmi ttances  in  the  present 
study. 

January  22,  1975 

The  area  of  the  United  States  east  of  the  Rocky  Mountains  was 
dominated  by  zonal  flow  with  a cutoff  low  over  New  Mexico  at  500  mb.  The 
surface  had  high  pressure  centers  over  western  Virginia  and  central  Iowa 
divided  by  a weak  frontal  system.  Richmond,  Virginia  37.53  N - 77.27  W 
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was  located  in  a general  area  of  scattered  fog  and  overcast  cirrus.  The 
location  had  an  observation  of  10/10  cirrus  from  4.4  - 5.7  km  and  10/10 
cirrostratus  from  6.7  - 8.7  km  and  5/10  cirrus  from  8.7  - 10.7  km.  The 
temperature,  pressure  and  water  vapor  profiles  were  obtained  from  the  Richmond 

1200Z  sounding.  The  cirrus  cloud  ice  content  was  6.4  gm  m-^  for  the  top 

7 -7 

layer,  67.8  gm  m for  the  middle  layer  and  104.2  gm  m for  the  lower 

layer  in  the  VTPR  field  of  view.  The  time  of  the  observation  and  satel- 
lite pass  time  was  1417Z.  The  SR  averages  were  248. 4°K  and  134.5  for  the 
IR  and  visible  channels,  respectively.  The  temperature  reading  is  the 
same  temperature  as  located  at  6.5  km  in  the  sounding. 

Fiqure  3.4  illustrates  comparisons  between  the  observed  and  com- 
puted upward  radiances  for  the  three  cirrus  cloud  cases  described  previously. 

Note  that  the  VTPR  scan  angles  are  denoted  in  Figure  3.3.  For  the  case 
when  thick,  layered  cirrus  clouds  are  present,  the  top  layer  consists  of 
only  50%  cloud  cover.  Thus,  calculations  were  made  by  averaging  the  re- 
sulting values  for  a three-  and  a two-layer  cloud  system.  The  calculated 
radiances  for  channels  4-8  are  consistently  lower  than  those  observed  from 
the  N0AA  4 satellite.  Possible  explanations  are  that  the  observed  value  of 
the  ice  content  was  overestimated,  or  the  particle  size  used  in  calculations 
was  too  small , or  the  cloud  temperature  was  improperly  assumed.  Comparisons 
between  the  observed  and  computed  radiances  for  the  single  cirrus  cloud  case 
reveal  systematically  higher  values  based  on  theoretical  simulations.  Again 
it  is  perhaps  caused  by  one  of  several  possibilities  such  as  an  underesti- 

* 1 

mation  of  the  vertical  ice  content  within  the  field  of  view  of  the  VTPR 
instrument.  In  our  judgement,  the  systematic  behavior  of  the  calculated 
and  observed  radiance  values  for  those  two  cases  indicates  to  certain  ex- 
tents the  reliability  of  the  theoretical  model  for  the  transfer  of  spectral 
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Figure  3.4.  Comparisons  between  satellite  observed  and  theoretically 

simulated  upward  radiances  for  the  three  cirrus  cloud  cases 
described  in  Figure  3.3. 
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infrared  radiation.  The  final  case  contains  an  overcast  lower  strato- 
cumulus  and  a light  cirrus.  It  is  assumed  that  the  low  cloud  is  a 
black  body  so  the  transfer  calculations  start  from  the  cloud  top.  Com- 
parisons with  observed  radiances  show  good  agreement  except  for  channel 
8.  In  order  to  find  the  physical  reason  for  the  deviation  in  this 
particular  channel,  careful  examinations  were  carried  out  for  transfer 
calculations.  The  observed  radiance  in  channel  8 is  about  98.3  watts 
cm"2  ster.  (cm'^)“^.  However,  utilizing  the  observed  temperature 
and  mixing  ratio  profiles  denoted  in  Figure  3.3a,  it  was  found  that 
the  clear  column  radiance  at  34.6°  zenith  angle  was  only  about  99.4 
watt  cm"2  ster.  (cnT^)-^.  In  view  of  the  considerable  reduction 
of  radiances  in  all  other  channels,  it  seems  that  the  observed  value  in 
channel  8 should  have  been  on  the  order  of  about  90  watt  cm  "2  ster. 
(crn*^)"1.  We  note  here  that  this  is  the  only  unsatisfactory  channel  com- 
parison. 

4> 

Within  the  uncertainties  involved  in  the  measured  radiances  and 
atmospheric  and  cloud  parameters,  and  the  assumptions  employed  in  the 
theoretical  model,  agreement  between  satellite  observed  and  theoreti- 
cally simulated  radiances  appears  satisfactory  for  these  three  selected 
cirrus  cloud  cases.  These  comparisons  establish  our  confidence  in  the 
theoretical  model  for  the  transfer  of  spectral  infrared  radiation.  It 
should  be  noted  that  since  the  peaks  of  the  weighting  functions  in  channel 
1-3  are  all  above  the  possible  cirrus  cloud  location,  cloud  effects  on  the 
upward  radiances  in  these  channels  can  be  ignored. 

3.3  Theoretical  Results  and  Cloud  Scene  Discussions 


3.3.1  Spectral  transmission  and  reflection  of  cirrus.  In  this  section 


we  present  some  results  of  the  radiative  properties  of  cirrus  clouds  and 
cirrus  cloudy  atmospheres . It  is  convenient  to  define  the  spectral  trans- 
mission Ta^  and  reflection  RA^j  of  a cloud  layer  as  follows  (see  Eqs.  (2.5) 
and  (2.6)): 


0*t-» ' Ov) 


(3.6) 


Physically,  the  transmission  defined  above  represents  the  actual  percent- 
age attenuation  of  upwelling  radiances  reaching  the  cloud  base,  and  it 
may  be  employed  for  the  parameterization  of  cloud  radiative  properties 
from  a satellite  point  of  view.  Note  that  both  transmission  and  reflec- 
tion include  contributions  of  cloud  emission,  which  is  automatically 
generated  in  transfer  calculations. 

Figure  3.5  illustrates  the  spectral  reflection  of  cirrus  as  a 
function  of  the  thickness  and  the  vertical  ice  content  for  VTPR  channels. 

The  left  and  right-hand  graphs  are  for  a tropical  (wet)  and  a midlatitude 
winter  (dry)  atmosphere,  respectively  (McClatchey  et  al . , 1971).  The  verti- 
cal scale  applies  to  the  lower-most  curves.  The  scales  for  other  curves 
are  to  be  obtained  by  subtracting  consecutively  a factor  of  0.2  such 
that  the  horizontal  bar  on  each  curve  is  zero.  Cloud  infrared  reflection 
appears  independent  of  the  atmospheric  profile.  Its  value  increases 
asymptotically  with  increasing  cloud  thickness.  Cloud  reflection  arises 
primarily  from  scattering  of  cloud  particles.  The  single-scattering  albedos 
are  about  the  same  for  channels  4,  5,  6,  and  8 resulting  in  similar  re- 
flection patterns.  For  channel  7,  however,  a higher  single-scattering 


albedo  is  derived  (see  Table  3.2).  This  results  in  larger  reflection  values 
as  evident  in  Figure  3.5.  Note  that  reflection  defined  earlier  includes 
emission  contribution  from  clouds,  primarily  the  warmer  cloud  base.  For 
a 1 km  thick  cirrus,  we  see  that  the  reflection  value  already  reaches 
about  0.4. 

Figure  3.6  shows  the  spectral  transmission  of  cirrus.  As  in 
Figure  3.5,  the  vertical  scale  applies  to  the  lower-most  curves.  The 
scales  for  other  curves  are  to  be  obtained  by  subtracting  successfully  a 
factor  of  0.2  such  that  the  horizontal  bar  on  each  curve  is  1.  Cloud 
infrared  transmission  decreases  with  increasing  thickness  and  reaches  a 
value  of  about  0.3  for  a cloud  thickness  of  5 km  for  channels  4-6  and  8. 

As  for  channel  7,  a cloud  whose  thickness  is  5 km  transmits  about  43%  of 
the  incident  upwelling  radiance  owing  to  a larger  single-scattering  albedo 
at  a wavenumber  of  535  cm"^.  Inspection  of  the  transmission  curves  indi- 
cates a similar  behavior  for  the  C02  and  window  channels.  In  particular,  ' 

\ 

variations  of  the  spectral  transmission  of  C02  channels  for  all  the  thick- 
nesses presented  here  are  less  than  about  3%.  Even  the  transmission  value 

i 

of  the  window  channel  does  not  vary  by  more  than  10%  with  respect  to  those 
of  the  C02  channels.  Larger  deviations  are  seen  to  take  place  for  larger 
thicknesses.  Since  cirrus  clouds  are  not  normally  thick,  we  may,  to  a - j 

good  approximation,  assume  that  the  cloud  transmission  values  are  con-  j 

stants  for  the  C02  and  window  channels. 

In  view  of  the  spectral  transmission  defined  in  Eq.  (3.6),  we 
may  parameterize  radiative  transfer  processes  in  cloud  layers  and  express 
the  upwelling  radiance  at  the  top  of  the  atmosphere  in  the  form 
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TROPICS  A k(  \ MIDLATITUDE  WINTER 


VERTICAL  ICE  CONTENT  (qm  cm'2  x I0'3) 

igure  3.6.  Spectral  transmission  of  cirrus  (see  text  for  definition)  as  a function  of  the  cloud 
thickness  for  a tropical  and  a mid-latitude  atmosphere.  The  lower  abscissa  denotes 
the  ice  content  scale  corresponding  to  the  cloud  thickness  employing  a particle 
number  density  of  0.05  cm"3. 


I.  = tB.(Ts)T.(zb,0)  + /0bB.[T(z)]d  T.(zb,z)} 


x TiC(2t,zb)Ti(-.,Zt)  + / B.[T(z)]d  T.(~,z). 


(3.7) 


Here  the  subscript  Av  is  replaced  by  i,  zb  and  zt  denote  the  cloud  base 
and  top  heights,  respectively,  and  T$  is  the  surface  temperature.  The 
first  term  represents  the  upwelling  radiance  reaching  the  cloud  base  and 
it  is  attenuated  by  the  cloud  and  gases  above  the  cloud  subsequently,  while 
the  second  term  denotes  the  absorption  and  emission  contributions  from  gases 
above  the  cloud.  According  to  spectral  transmission  calculations  reported 
previously,  T-jC  (i  = 4,5,6  and  8)  may  be  considered  as  one  unknown  para- 
meter which  is  composed  of  the  cloud  ice  content  (or  optical  depth)  infor- 
mation. Thus,  for  overcast  cirrus  cloud  cases,  it  seems  feasible  to  con- 
struct numerical  procedures  for  recovering  the  cloud  transmission  (in  terms 
of  ice  content,  for  example)  and,  perhaps,  temperature  profile  simultaneously 
utilizing  CO2  and  window  channels. 

3.3.2  Upward  radiances  and  possible  ice  content  determination.  The  most 
important  quantity  associated  with  satellite  sensing  of  the  atmosphere  is 
the  upward  radiance  at  the  top  of  the  atmosphere.  Upward  radiances  are  the 
only  information  for  the  inference  of  the  structure  and  composition  of  the 
atmosphere. 

Figure  3.7  shows  graphs  of  upward  radiances  in  VTPR  channels  as 
functions  of  the  zenith  angle  for  various  cloud  ice  contents  in  a tropical 
atmosphere.  For  a given  ice  content  value,  it  is  clear  that  upward  radiances 
increase  from  channel  4 to  channel  8 (note  that  upward  radiances  in  channel 
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ZENITH  ANGLE  (dag) 

Figure  3.7.  Upward  radiances  at  the  top  of  the  atmosphere  in  VTPR  CO2  and  window 
channels  as  functions  of  the  zenith  angle  for  various  cloud  ice  con- 
tents ranging  from  14  to  112  gm  nr  2.  The  model  atmosphere  used  in 
calculations  is  a tropical  profile. 

40 


7 is  not  shown).  This  is  probably  due  to  decreasing  absorption  strength 


from  the  C02  band  center  to  the  wing.  For  an  ice  content  of  112  gm  m , 

! we  see  that  upward  radiances  in  channels  4-6  are  about  the  same  with  pat- 

i terns  close  to  isotropic  radiation  except  in  the  limb.  Upward  radiances 

in  these  channels  show  strong  limb  darkening  effect  for  zenith  angles 
greater  than  about  30°  or  so.  From  0°  (nadir)  to  about  30°  zenith  angles, 
upward  radiances  do  not  vary  appreciably.  With  respect  to  each  channel, 
increasing  the  ice  content  of  cirrus  leads  to  a reduction  of  upward  radi- 
ances. This  is  particularly  evident  for  channel  8.  Physically,  this 
means  that  clouds  affect  more  significantly  radiance  measurements  in  the 
< window  channel  as  compared  with  the  CO2  channels.  Spectral  transfer 

calculations  were  also  carried  out  for  other  climatological  profiles 
tabulated  by  McClatchey  et  al . (1971).  These  calculations  reveal  that 
relationships  between  upward  radiances  and  ice  contents  depend  somewhat 
on  the  atmospheric  temperature  profile,  except  for  very  large  values  of 
ice  contents  which  dominate  transfer  processes. 

In  Figure  3.8,  we  present  the  upwelling  (nadir)  radiance  of 
VTPR  channels  for  two  model  cirrus  cloudy  atmospheres  representing  wet 
(tropics)  and  dry  (midlatitude  winter)  conditions,  respectively.  The 
solid  curves  in  this  figure  denote  upwelling  radiances  for  clear  columns. 
For  the  tropical  atmosphere,  systematic  reductions  of  the  upwelling  radi- 
ance due  to  the  increase  of  the  ice  content  are  i 1 lustrated.  In  reference 
to  Figure  3.1,  it  is  seen  that  the  peak  of  the  weighting  function  for 
channel  4 is  located  at  about  the  position  of  cirrus  clouds.  Consequently, 
even  a thin  cirrus  produces  a significant  signature  of  radiance  reduction. 
In  the  water  vapor  rotational  band  (smaller  wavenumber)  larger  upwelling 
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VTPR  CHANNEL 

Figure  3.3.  Upwelling  (nadir)  radiances  of  VTPR  channels  for  a trooical  (left-hand  side)  and  a mid- 
latitude  (right-hand  side)  atmosphere.  The  solid  curves  denote  upwelling  radiances  for 
clear  columns.  Systematic  reductions  caused  by  the  increase  of  the  cirrus  ice  content 
is  shown  by  dashed  curves. 


radiances  arise  from  the  behavior  of  the  Planck  function  and  the  emission 
from  water  vapor.  The  systematic  reduction  caused  by  the  increase  of  ice 
content  is  so  distinct  that  the  observed  radiances  obtained  from  VTPR 
channels  in  a completely  cloudy  atmosphere  may  be  utilized  to  evaluate 
the  amount  of  ice  within  the  f ield-of-view.  This,  of  course,  would  rely 
upon  a first  estimation  of  clear  column  radiances  from  available  atmospheric 
profiles.  Upwelling  radiances  in  a midlatitude  winter  atmosphere  exhibit 
similar  reduction  patterns,  but  their  values  are  lower  owing  to  the  lower 
temperatures  in  the  troposphere.  Reauction  values  in  cases  when  ice  con- 

O 

tents  are  84  and  112  gm  m are  very  close,  reveal ing  that  the  cloud  is 
reaching  a black-body  thermodynamic  state. 

Having  the  upward  radiance  calculated  for  a series  of  cirrus  ice 
content  values  in  various  model  atmospheric  profiles,  the  following  proce- 
dures seem  feasible  for  the  ice  content  determination.  Assume  that  within 
the  field-of-view  of  VTPR  channels  there  is  n portion  of  cirrus  cloudiness, 
the  upward  radiance  in  a partly  cloudy  atmosphere  is  given  by  Eq.  (2.18) 
in  the  form  (av  -*•  i ) . 

IiPC  = n 1^  + (l-n)I-j^C.  (3.8) 

It  can  be  shown  that  for  any  two  cloud-contaminated  radiances  at  distinct 
wavenumbers  that 

j PC  _ i NC  / j C _ j NC 

/ PC  \ NC  / j"1  " C 1 , NC  = (3,9) 

7+1  “ 7 + 1 / 7+1  " 7 + 1 
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where 

T C t PC  T NC  4 _ 4 - c 7 o 
I ^ ^ I ^ ^ 1^  , 1 bj  6(  /)  o. 

As  mentioned  previously,  clear  column  radiances  are  to  be 
estimated.  Estimation  may  be  accomplished  by  employing  available  atmos- 
pheric profile  analyses  such  as  those  provided  by  NMC  or  from  climatology. 

Although  Eq.  (3.9)  now  contains  two  unknown  parameters  1^  and  these 

two  variables  are  both  functions  of  the  ice  content.  Based  on  transfer 

l 

calculations,  for  a given  ice  content  in  a specified  atmosDhere,  both 
values  may  be  derived.  Upon  inserting  a series  of  ice  content  values, 
we  may  obtain  a best  estimate  such  that  Eq.  (3.9)  holds.  This  procedure 
can  be  carried  out  for  pairs  of  adjacent  channels  in  which  the  optical 
property  of  ice  varies  insignificantly.  Consequently,  each  pair  gives  an 
ice  content  estimation.  By  further  inspection  of  these  estimations,  a final 
value  for  ice  content  which  satisfies  all  the  channel  observations  could  be 

1 

determined.  This  could  be  accomplished  by  examining  the  mean  and  standard 
deviation  with  respect  to  1 according  to  Eq.  (3.9)  for  all  available  channel 
observations.  Once  an  ice  content  value  has  been  chosen,  I-jC  is  known  and 
from  Eq.  (3.8),  the  cloud  cover  information  can  be  calculated.  An  averaged 
value  of  n may  be  subsequently  obtained  such  that  it  satisfies  all  the  channel 
observations. 

The  above  procedures  were  followed  for  the  three  cirrus  cloud  cases 
described  in  section  3.2.  The  best  ice  content  estimates  that  we  obtained 
were  10,  8,  and  112  gm  m" 2 for  11  January,  1974,  1 February,  1974  and  22 
January,  1975,  respectively,  whereas  the  in-si tu  aircraft  observations 
gave  11,  2 and  176  gm  m-S  respectively.  Since  these  are  the  only  three 
available  cirrus  cloud  cases,  it  is  not  possible  to  perform  significant 
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statistical  error  analyses. 

It  should  be  emphasized  that  the  semi -empi rical  method  for  ice 
content  estimation  makes  use  of  the  results  from  transfer  calculations  and 
assumes  a prior  knowledge  of  the  clear  column  radiance.  The  latter  assump- 
tion relies  upon  reliable  atmospheric  profiles  which  are  unknown  in  cloudy 
conditions.  Hence,  the  best  estimation  for  the  profile,  which  is  needed, 
may  lead  to  the  uncertainty  in  the  ice  content  determination.  Nevertheless, 
the  procedure  described  here,  at  least,  represents  an  objective  and  workable 
means  for  recovering  cloud  information  from  satellite  sensing. 

3.4  Conclusions 

The  infrared  transfer  model  described  in  Section  2 was  applied 
to  VTPR  channels  of  the  NOAA  4 satellite  to  simulate  upward  radiance  at 
the  top  of  the  atmosphere  in  cirrus  cloud  conditions.  Available  transmit- 
tances  for  each  channel  were  employed  and  modified  to  incorporate  them 
into  the  spectral  infrared  transfer  program.  Comparisons  between  satellite 
observed  and  theoretically  simulated  radiances  reveal  systematic  agreement 
for  selected  cirrus  cloud  cases  involving  thin  and  thick  cirrus,  and  thin 
cirrus  with  stratocumulus  below.  Theoretical  calculations  employ  observed 
atmospheric  profiles  from  radiosondes,  and  the  best  estimated  cloud  para- 
meters from  aircraft  observations  under  the  satellite  pass.  These  compari- 
sons give  confidence  in  the  theoretical  model  for  the  spectral  infrared 
transfer  in  cloudy  atmospheres. 

Utilizing  the  spectral  infrared  transfer  model,  radiative  proper- 
ties of  cirrus  cloudy  atmospheres  in  VTPR  channels  were  investigated  for 
two  model  atmospheric  profiles.  The  spectral  transmission  and  re- 
flection, defined  as  the  ratios  of  the  upwelling  radiance  at  the  cloud 
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bottom  to  that  at  the  cloud  top  and  to  the  downwelling  radiance  at  the 
cloud  bottom,  respectively,  were  calculated.  We  show  that  effect  of  the 
atmospheric  profile  on  cloud  transmission  and  reflection  is  insignificant 
and  that  close  resemblance  of  these  two  parameters  is  found  for  VTPR 
channels  4-6  and  8.  We  further  demonstrate  that  using  the  concept  of 
cloud  transmission  it  appears  possible  to  parameterize  complicated  radi- 
ative transfer  processes  through  clouds. 

Calculations  of  upward  radiances  were  made  for  several  model 
cirrus  cloudy  atmospheres  by  varying  the  value  of  cloud  ice  content.  For 
a given  atmosphere,  it  is  shown  that  the  decrease  of  upward  radiances  is 
closely  related  to  the  ice  content  of  cirrus  clouds.  Systematic  reduc- 
tion patterns  of  the  upwelling  radiance  caused  by  the  increase  of  the 
ice  content  are  clearly  illustrated  for  VTPR  channels  employing  tropical 
and  midlatitude  atmospheric  profiles.  Finally,  procedures  are  outlined 
for  estimating  the  cirrus  ice  content  in  partly  covered  conditions.  Em- 
ploying three  cirrus  cloud  cases  in  which  satellite  observed  radiances 
are  available,  we  demonstrate  that  the  ice  content  and  the  amount  of 
cloud  cover  may  be  estimated  objectively. 


SECTION  4 


APPLICATION  TO  NIMBUS  VI  HIRS  CHANNELS 

4.1  Characteristics  of  HIRS  Channels  and  Model  Atmosphere 

This  section  presents  a method  for  cloud  type  and  cloud  mass 
determination  utilizing  the  High  Resolution  Infrared  Sounder  (HIRS)  data 
based  upon  parameterization  of  calculations  from  the  spectral  infrared 
model  described  in  Section  2.  The  Nimbus  VI  HIRS  instrument  is  a third 
generation  infrared  radiation  sounder  and  is  similar  to  the  Infrared  Tem- 
perature Profile  Radiometer  (ITPR)  on  the  Nimbus  V Satellite.  The  instru- 
ment scans  perpendicular  to  the  satellite  subtrack.  There  are  42  scan 
spots  per  scan  line  with  a resolution  of  23  km  near  nadir  and  31  km  at 
the  extremes  of  the  scan.  The  Nimbus  VI  satellite  was  successfully  launched 
in  June  1975.  The  HIRS  instrument  had  short  periods  during  which  all 
channels  were  operating  successfully.  We  have  obtained  a sample  of  good 
data  from  Smith  (private  communication).  The  data  correspond  to  a period 
of  20-30  August,  1975  and  cover  a geographical  area  from  80°-150°W  and 
20°-50°N.  The  HIRS  instrument  senses  infrared  radiation  in  17  channels 
which  include  7 channels  in  the  15  pm  CO2  band,  5 channels  in  the  4.3  um 
CO2  band,  water  vapor  channels  at  6.8  pm,  and  3 channels  in  windows  at 
11  pm,  3.68  pm  and  0.69  pm. 

The  data  we  received  were  not  corrected  for  solar  radiation  con- 
tamination in  the  short  wave  channels.  Reflected  sunlight  has  a pronounced 
affect  on  channel  16,  the  short  wave  window,  which,  however,  was  not  used 
in  the  present  analysis. 


The  information  compiled  in  Table  4.1  and  Figure  4.1  are  taken 
from  Smith,  et  al.  (1975).  Figure  4.1  depicts  weighting  functions  for  the 
channels  in  Table  4.1.  The  peak  in  this  figure  indicates  the  approximate 
location  in  the  troposphere  from  which  its  energy  is  derived.  Only  those 
channels  whose  weighting  functions  peak  below  100  mb  would  be  substantially 
affected  by  clouds. 


Table  4.1.  HIRS  channel  characteristics. 
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Figure  4.1 

The  weighting  functions  of  the 

HIRS  channel 

The  model  atmosphere  used  in  the  theoretical  analysis  was  the 
mid-latitude  summer  atmosphere  described  by  McClatchey  et  al.,  (1972). 

Cloud  types  were  restricted  to  middle  (al tostratus)  and/or  high  (cirrus). 

The  model  atmosphere  was  divided  in  such  a way  that  it  would  coincide  with 
the  40  pressure  levels  used  in  the  clear  column  radiance  program  (CCR) 
developed  at  NOAA/NESS,  and  supplied  to  us  for  our  use.  The  program 
utilizes  predetermined  transmission  profiles  which  can  be  adjusted  as  a 
function  of  the  atmospheric  temperature  profile  used.  It  includes  contri- 
butions from  all  the  major  gases  listed  in  Table  4.1.  To  further  utilize 
this  predetermined  vertical  structure  and  facilitate  the  execution  of  the 
cloud  transfer  program,  the  base  of  the  middle  cloud  was  fixed  at  700  mb  and 
the  top  of  the  high  cloud  was  fixed  at  250  mb.  The  thicknesses  of  the  high 

and  middle  clouds  were  allowed  to  vary  from  1.26  to  4.12  km  and  from  0.35 

to  2.68  km,  respectively  in  the  theoretical  sensitivity  analyses.  Thus, 
two  distinct,  realistic  cloud  layers,  with  changeable  thicknesses,  are 
formed  and  the  model  is  developed  to  utilize  all  the  possible  thickness 
variations. 

The  single  scattering  parameters  for  altostratus  and  cirrus  cloud 
were  calculated  by  using  scattering  computations  for  the  central  wavenumber 
in  each  band.  The  cirrus  cloud  was  assumed  to  be  entirely  ice  and  the  alto- 
stratus  cloud  was  assumed  to  be  all  water.  Table  4.2  lists  the  optical  pro- 
perties of  the  ice  cylinders.  The  real  and  imaginary  parts  of  the  index  of 
refraction  are  taken  from  the  recent  measurement  by  Schaaf  and  Williams  (1973). 

From  Table  4.2  it  is  apnarent  that  the  real  part  of  the  index  of  refraction 

increases  with  decreasing  wavenumber  and  the  imaginary  part  decreases  with 
increasing  wavenumber.  The  single  scattering  albedo  is  largest  in  the 


window  channels  (channels  16  and  8)  and  the  extinction  cross  section  is 
greatest  in  the  15  pm  band  (channel  4-7).  Table  4.2  also  contains  the 
optical  properties  of  the  water  droplets.  The  real  and  imaginary  parts 
of  the  index  of  refraction  for  each  band  were  taken  from  Hale  and  Querry 
(1973).  From  Table  4.2  it  can  be  seen  that  the  real  part  of  the  index 
of  refraction  generally  increases  with  increasing  wavenumber  and  the 
imaginary  part  decreases.  The  single  scattering  albedo  is  the  largest 
at  the  largest  wavenumber  and  the  extinction  cross  section  is  the  greatest 
in  the  water  vapor  bands  (channels  9 and  10).  Only  one  calculation  was 


made  for  each  of 

the  C02  bands 

since 

the  real  and  imaginary  parts 

of  the 

index  of 

refraction  vary  slightly. 

Table  4 

.2.  Optical  properties  of  cirrus  and  altostratus 

clouds. 
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.2  Ratioing  of  Upwelling  Radiances:  A Case  Study  on  the  Cloud 

Thickness  Estimation 

For  each  channel,  the  theoretical  model  calculated  upwellinq 
clear  column  and  cloudy  radiances  at  the  different  scan  angles  of  the  HI RS 
instrument,  given  an  atmosphere,  cloud  temperature  and  cloud  composition 
for  each  cloud  type.  In  addition,  for  computational  purposes  the  cloud 
top  height  for  the  cirrus  and  the  cloud  base  for  the  middle  cloud  were 
held  constant.  In  the  real  atmosphere,  the  temperature  profile  is  always 
changing,  the  cloud  composition  is  variable,  and  the  location  of  the  clouds 
in  the  vertical  is  never  constant.  A method  that  will  minimize  the  effects 
of  the  model  assumptions  when  parameterizing  the  real  atmosphere  is  desir- 
able. To  minimize  the  computer  requirements  and  redundancy  in  analysis 
all  the  subsequent  relationships  were  derived  for  a scan  angle  of  0°. 

These  same  relationships  could  be  developed  for  any  scan  angle  or  range 
of  scan  anges. 

To  minimize  the  model  assumptions  and  to  normalize  the  real  data 
for  the  model,  the  clear  column  radiances  were  divided  into  the  cloudy  ra- 
diances for  each  channel.  This  ratio  represents  the  relative  reduction  of 
upwelling  radiances  due  to  the  cloud  effects  in  the  atmosphere  and  the 
quantity  is  dimensionless.  There  are  several  advantages  in  working  with 
these  ratios.  Physically,  the  ratio  will  reduce  the  effect  of  the  change 
in  atmospheric  profile  on  the  cloudy  radiances.  In  addition,  the  degrada- 
tion of  the  ratio  coupled  with  the  peaking  of  the  weighting  function  gives 
an  immediate  indication  of  clouds  at  that  level  or  above.  Another  advant- 
age of  the  ratio  technique  is  that  in  addition  to  normalization  between 
different  atmospheres,  the  effects  of  clouds  on  the  ratio  of  each  channel 
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can  be  compared  with  other  channels.  In  the  subsequent  analysis,  channels 
4-7  of  15  uni  CO-,  band,  11-14  of  the  4.3  um  CO2  band,  8 at  11.11  um  in  the 
window,  9 in  the  8.2  um  water  vapor  band  and  10  in  the  6.3  um  water  vapor 
band  were  utilized. 

The  theoretical  model  was  executed  for  cirrus  thicknesses  of  1, 

2,  3,  and  4 km  for  the  case  study  described  below.  These  ratios  are  plotted 
in  Figure  4.2  (solid  lines).  Examination  of  Figure  4.2  shows  that  cirrus 
is  transparent  in  the  infrared  spectrum.  There  is  a successive  decrease  of 
upwelling  radiances  in  all  channels  as  cloud  thicknesses  increase.  The 
successive  decrease  of  the  ratios  continues  all  the  way  to  4 km.  Between 
3 and  4 km  the  decrease  in  ratio  in  each  channel  is  only  slight  indicating 
that  near  4 km  the  cirrus  cloud  used  in  the  model  is  becoming  opaque  to  the 
infrared  radiation. 

Comparison  of  the  channels  in  Figure  4.2  gives  a good  indication 
of  the  effects  of  clouds  on  upwelling  radiances  in  different  wavenumber 
regions  of  the  infrared  spectrum  relative  to  the  peak  of  the  channel's 
weighting  function.  The  ohysical  factors  that  influence  the  ratio  for  a 
given  cloud  type,  cloud  thickness  and  channel  could  be  determined  from 
ratios  presented  in  Figure  4.2.  In  the  present  model,  the  cloud  top  is 
constant  at  250  mb.  Using  a 2 km  cirrus  as  an  example,  the  interaction  of 
these  two  physical  factors  (the  weighting  function  and  the  channel  wavenumber) 
are  examined.  The  channels  that  are  in  the  15  um  CO2  band  show  a decreasing 
ratio  as  the  weighting  functions  peak  deeper  in  the  atmosphere.  The  same 
decrease  is  noted  in  4.3  um  CO2  channels  but  to  a greater  degree.  This  fact 
would  indicate  a greater  impact  of  clouds  on  shorter  wavelengths.  For 
example,  channel  7 (peak  at  900  mb)  has  a ratio  nearly  double  that  of 
channel  12  (peak  at  850  mb). 
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Figure  4.2.  Comparison  of  theoretical  results  (solid  curves)  and  satellite  HIRS  data  □ , o >•) 

for  the  11  HIRS  channels  in  terms  of  the  cloudy  to  clear  column  radiance  ratios. 


There  were  26  passes  available  in  the  20-30  August,  1975,  HIRS 
data  set.  Inspection  of  NOAA  4 mapped,  normalized,  and  qridded  mosaics 
from  the  Very  High  Resol ution  Radiometer  (VHRR)  for  this  time  period  re- 
veals one  pass  of  particular  interest.  On  25  August  a cirrus  shield  was 
associated  with  a squall  line  near  40°N  in  the  central  U.S.  The  corres- 
ponding HIRS  pass  was  analyzed  for  cloud  information  at  scan  angles  closest 
to  nadir.  Figure  4.3  shows  both  the  broad  band  visible  channel  (0.5-0. 7 um) 
and  the  broad  band  window  channel  (10.5-12.5  pm). 

To  apply  the  theoretical  results  to  the  real  HIRS  data  a method 
had  to  be  developed  to  reduce  the  cloudy  radiances  to  ratios.  The  theo- 
retical and  the  real  data  are  assumed  compatable  since  both  atmospheres 
are  mid-latitude  summer.  To  obtain  a clear  column  radiance  that  is  re- 
presentative of  a localized  area,  the  HIRS  data  was  examined  in  conjunction 
with  the  satellite  pictures.  The  point  chosen  for  the  clear  column 
radiance  (at  about  38°N,  91 °W)  was  then  used  to  obtain  the  cloudy  ratios 
for  that  day.  The  point  chosen  has  the  same  scan  angle  as  the  cloudy 
radiances. 

In  reference  to  Figure  4.3,  we  have  selected  a number  of  loca- 
tions at  about  40°N  and  90°W  along  the  satellite  track  for  the  inference 
of  cloud  thickness  based  on  the  theoretical  sensitivity  analyses  described 
previously.  Examinations  of  these  cloud  photographs  as  well  as  of  synoptic 
surface  and  500  mb  maps  show  that  the  cloud  is  cirrus  with  no  low  cloud 
present  below.  Six  scan  radiance  ratios  are  depicted  in  Figure  4.2.  The 
scan  points  at  38. 7N/91 .2W(v)  and  38. 7N/91 . 6W( 0 ) represent  areas  that  are 
close  to  the  cirrus  edge.  Based  on  the  comparisons  with  the  theoretical 
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curves,  thicknesses  of  0.4  and  1 km  were  estimated.  As  the  satellite  scan 
goes  deeper  into  the  cloud  center,  a thickness  of  greater  than  4 km 
( 40 . 1 N/92 . 1 W , A ) is  obtained.  Other  scans  further  north  under  investiga- 
tion indicate  cirrus  thicknesses  on  the  order  of  2 to  3 km  as  illustrated 
by  square,  hexagon  and  black  dot.  Note  that  the  resolution  of  HIRS  in- 
struments is  approximately  23  km  at  nadir.  From  this  case  study  it  appears 
that  the  HIRS  radiance  ratios  resemble  reasonably  well  those  calculated 
from  the  theoretical  analyses.  The  4.3  and  15  vm  C02  channel  behave 
especially  well.  However,  the  water  vapor  channels  show  some  non-system- 
atic  variation,  perhaps  owing  to  the  use  of  the  observed  clear  column 
radiance  which  arises  from  an  atmosphere  differing  from  cloudy  cases. 

4.3  Parameterization  of  Radiance  Calculations  from  Cloudy  Atmospheres 

4.3.1  Cloud  type  determination.  To  utilize  the  theoretical  results  in 

a real  atmosphere,  we  first  need  to  determine  the  vertical  location  of  the 
cloud  in  the  atmosphere.  This,  in  effect,  will  give  a measure  of  the  cloud 
type.  The  theoretical  calculations  were  carried  out  for  cirrus  thicknesses 
of  1,  2,  3,  and  4 km  and  middle  cloud  thicknesses  of  0.3,  0.5,  0.7,  1.0, 
2.0,  and  3.0  km.  Examination  of  the  upwelling  ratios  shows  that  the 
transparent  quality  of  cirrus  is  evident  in  the  infrared  spectrum.  There 
is  a successive  decrease  of  upwelling  radiances  in  all  channels  as  the 
cloud  thicknesses  increase.  The  successive  decrease  of  the  ratios  con- 
tinues all  the  way  to  4 km.  Between  3 and  4 km  the  decrease  in  ratio 
in  each  channel  is  only  slight  indicating  that  near  4 km  the  cirrus  cloud 
used  in  the  model  is  becoming  opaque  to  the  infrared  radiation,  and  the 
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cloud  temperature  now  becomes  important  in  the  transfer  process.  At 
cirrus  thicknesses  greater  than  4 km,  the  upwelling  radiance  is  represent- 
ative to  the  Planckian  temperature  of  the  cloud  top.  Analysis  of  middle 
cloud  thickness  indicates  the  cloud  mass  becomes  opaque  between  2 and  3 km. 
This  may  be  attributed  to  greater  particle  concentrations  found  in  middle 
clouds  compared  to  cirrus  clouds. 

For  the  purpose  of  determining  the  cloud  type  from  the  theo- 
retical data,  the  channels  were  rearranged  in  order  of  decreasing  ratios 
for  a moderately  thick  cirrus.  The  rearrangement  of  the  channels  is 
shown  in  Figure  4.4.  Analysis  of  this  figure  shows  the  effects  of  wave- 
number  and  weighting  function  peak.  The  long  wave  channels  4-7  all  have 
greater  ratios  than  the  short  wave  channels  11-14,  even  though  they  have 
similar  weighting  function  peaks,  which  would  indicate  that  the  cloud 
mass  has  more  impact  on  shorter  wavelengths.  The  location  of  channel 
14  between  channels  5 and  6 would  indicate  that  there  is  an  overlap  be- 
tween  the  two  CO2  bands  depending  on  the  weighting  function  peak.  The 
locations  of  channels  10  and  9 are  according  to  weighting  function  peaks, 
with  channel  10  having  a higher  ratio  and  a higher  weighting  function 
peak  than  channel  9.  The  location  of  their  ratios  within  the  C0?  channel 
ratios  in  the  figure  is  due  to  the  model  atmosphere  humidity  profile. 

By  increasing  the  water  vapor  amount,  there  would  be  an  increase  in 
attenuation  and  thus  a lower  ratio.  Conversely,  a decrease  in  the  water 
vapor  would  decrease  the  attenuation  and  increase  the  ratio. 

The  radiance  ratios  and  channel  numbers  shown  in  Figure  4.4 
were  fitted  with  a linear  equation  by  means  of  regression  analysis. 
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Figure  4.4  The  best  straight  line  fit  for  four  cirrus  cloud 
tnicknesses  (Part  a)  and  for  six  middle  cloud 
thicknesses  (Part  b) . 
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Part  a is  for  cirrus  clouds  and  Part  b is  for  middle  clouds.  Each 
straight  line  is  labeled  with  the  thickness  it  represents.  The  equa- 
tion is  for  the  form  y = a-|R  + a0  where  y equals  the  number  assigned 
to  the  channel  (see  Figure  4.4)  and  R is  the  ratio  for  any  given  channel 
and  thickness.  These  parameters  are  defined  by 


ERy  - rREy/n 

= J U • 

ER^  - (eR)  /n 

(4.1) 

J j 

a0  = y - a-|  R, 

(4.2) 

y = sy/n, 

(4.3) 

j 

R = ER/n, 

(4.4) 

j 

and  the  square  of  the  correlation  coefficient  of  the  linear  regression 
analysis  is  given  by  (see  e.g.,  Panofsky  and  Brier,  1968) 

[ERy  - EREy/n]2 

r2  = -J-5 J-  J-2- o 2 ' (4.5) 

[ER  - (eR)  /n]  [Ey2  - (Ey) 2/n] 

j j j j 

In  the  equations  the  summations  are  for  all  11  channels. 

Table  4.3  contains  the  slope,  y-intercept,  and  r2  for  each  thick- 
ness of  cirrus  and  middle  clouds.  Analysis  of  the  table  shows  that  the  cor- 
relation (r2)  is  much  better  in  the  cirrus  case  (0.96)  when  compared  to  the 
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middle  cloud  case  (0.90).  The  decrease  in  correlation  in  the  middle  clouds 
and  the  1 km  cirrus  cases  can  be  explained  by  the  placement  of  the  ratio  of 
channel  10.  The  original  arrangement  was  done  for  a thick  cirrus  case 
where  the  effects  of  the  cirrus  cloud  mass  on  channel  10  were  present.  In 
the  middle  cloud  and  1 km  cirrus  cases  there  is  no  effect  of  the  cloud  mass 
on  channel  10.  Therefore,  the  location  of  channel  10  in  the  fitting  process 
decreases  the  correlation  for  the  1 km  cirrus  and  the  middle  cloud  cases. 

The  linear  fit  developed  here  will  be  used  when  parameterizing  actual  HIRS 
data  in  the  next  section.  Another  important  thing  to  note  in  this  analysis 
is  that  the  magnitude  of  the  slope  is  greater  than  the  magnitude  of  the  y- 
intercept  for  all  cirrus  thicknesses  and  the  reverse  is  true  for  middle 
clouds.  This  fact  will  be  used  in  subsequent  analysis  to  distinguish  be- 
tween middle  and  cirrus  clouds. 

Table  4.3.  Slope,  y-intercept,  and  correlation  coefficient 
squared  for  several  cloud  thickness  for  cirrus 
(part  a)  and  middle  (part  b)  clouds. 

Part  a (cirrus  clouds) 


1 

Thickness  (km) 

2 3 

4 

slope 

y-intercept 

-21 . 75 
21.69 

-16.82  -14.79 

15.70  13.68 

-13.94 

12.84 

r^ 

.872 

.968  .977 

.977 

Part  b (middle  clouds) 

.3 

Thickness  (km) 

.5  .7  1. 

2. 

3. 

slope 

y-intercept 

-22.78 

23.56 

-19.80  -18.80  -18.10 

20.70  19.70  18.80 

-17.50 

18.10 

-17.50 

18.10 

r2 

.904 

.925  .925  .913 

.906 

.906 
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4.3.2  Ice  and  water  content  determination.  A measure  of  liquid  water 
and  ice  content  distribution  over  the  planet  is  important  to  many  areas  of 
meteorology.  These  include  the  initialization  of  numerical  models  to  in- 
crease long  term  forecast  accuracy  and  the  study  of  global  circulation  and 
climatic  change  models.  From  the  theoretical  calculations,  each  cloud 
thickness  is  representative  of  a cloud  mass  since  only  one  drop  size  dis- 
tribution was  used  for  water  clouds  and  only  one  particle  concentration 

was  used  for  ice  clouds.  For  water  clouds  a liquid  water  content  of 

-3  -3 

0.15  gm  m was  assumed,  and  for  ice  clouds  an  ice  content  of  0.0283  q m m 

was  used. 

The  recovery  of  liquid  water  and  ice  content  from  the  fitted 

ratios  is  possible  by  knowing  within  certain  limits  the  theoretical  mass 

of  the  cloud  for  which  the  ratios  are  being  calculated.  In  the  cirrus  cloud 

_o 

case,  it  was  found  that  for  a cloud  mass  of  less  than  28.3  gm  m the  cor- 
relation coefficient  decreased  significantly  below  95  percent  and  the 

ratios  could  no  longer  be  fit  with  a straight  line.  This  was  also  true 

_2 

for  middle  clouds  with  a mass  of  less  than  45  gm  m . Thick  middle  clouds 
(greater  than  390  gm  m“2)  had  the  same  ratios  as  the  300  gm  m-2  (2  km 
thickness)  cloud.  This  would  indicate  that  the  cloud  mass  is  a black 
body  and  all  the  energy  from  below  is  absorbed  by  the  cloud  layer  and 
reradiated.  The  113.2  gm  m-2  (4  km  thick)  cirrus  still  seems  to  be 
somewhat  transparent,  but  this  cloud  mass  was  near  the  maximum  that  could 
be  handled  by  the  theoretical  model. 

To  recover  ice  or  water  content  from  the  parameterizations  of 
the  theoretical  radiance  ratios,  the  slopes  and  y-intercepts  for  the  two 
cloud  types  could  be  fit  by  a least  square  loqarithmic  function  that  is 
given  for  the  slope  by 
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al  = C1  + c2  &z 


(4.6) 


and  for  the  y-intercept  by 


aQ  = c3  + c4  In  az. 


(4.7) 


where  A z is  the  thickness  of  the  ice  or  water  cloud,  and  c^ , C£,  c-j,  and 
c4  are  constants  to  be  determined  for  each  set  of  slopes  and  y-intercepts. 
For  the  slope  ap  we  have 


C1  = 


£(a,ln&z)  - ElnAzEa,/n 

J j j 

E(lnAz)2  - (ElnAz)2/n 


(4.8) 


c2  = ^ (sa1  - cjElnAz). 
i i 


(4.9) 


These  equations  are  derived  by  substitution  of  the  natural  logarithm  of 
the  thickness  for  R, and  a1  for  y in  Eqs.  (4.1)  and  (4.2).  The  square  of 
the  correlation  coefficient  is  given  by 


[za^lnAz  - ElnAz  Ea^/n]* 


[E(lnAz)2  - (zlnAz)2/n]  [sa  2 - (Ea.)2/n] 
i i i 1 i 1 


(4.10) 


which  is  again  a direct  substitution  into  Eq.  (4.6).  The  y-intercept,  c3, 
c4  and  rc  can  be  calculated  by  substituting  a0  for  a1  in  Eqs.  (4.8)  - 
(4.10).  The  summation  is  over  six  cloud  thicknesses  for  middle  clouds 
and  over  four  cloud  thicknesses  for  cirrus. 

Using  the  data  from  Table  4.3,  the  slopes  and  y-intercepts  were 
fitted  for  each  cloud  type.  The  resulting  fits  are  shown  in  Figure  4.5 
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for  middle  clouds  and  Figure  4.6  for  cirrus  cloud.  For  middle  clouds  c-j 
and  c^  are  equal  to  -18.53  and  3.007,  respectively,  with  the  correlation 
coefficient  of  0.984,  while  c3  and  c^  equal  19.19  and  -4.089  with  a cor- 
relation coefficient  of  0.960.  For  the  cirrus  case  the  values  of  Cp  c 2> 
c3,  and  c^  are  -21.40,  5.75,  21.19,  and  -6.526,  respectively.  The  constants 
Cj  and  c^  and  a correlation  coefficient  of  0.988  and  C3  and  c4  had  a cor- 
relation of  0.983. 

Upon  substituting  a1  and  a2  in  Eqs.  (4.6)  and  (4.7)  into  Eq. 

(4.2),  the  thickness  of  any  middle  or  high  cloud  is  given  by 


az  = exp  [(y-Tfcj  - c3)/(c2T*  + c4)],  (4.11) 

The  cloud  mass  is  then  derived  by  multiplying  the  cloud  thickness  expressed 
in  km  by  28.3  gm  m' 2 for  cirrus  cloud  and  1.5  x 102  gm  m~2  km-1  for  middle 
clouds. 

4.4  Data  Description  and  Selection 

The  method  to  determine  cloud  type  and  ice  or  water  content 
developed  in  the  previous  section  was  applied  to  a selected  set  of 
Nimbus  VI  High  Resolution  Infrared  Sounder  (HIRS)  data.  The  results 
of  this  appl ication  were  compared  to  the  Air  Force  Three  Dimensional 
Nephanalysis  (3DNEPH)  and  satellite  cloud  pictures  since  these  were 
the  only  source  of  consistent  cloud  information  available  on  a routine 
basis. 

4.4.1  The  three-dimensional  nephanalysis  data.  The  Three-Dimensional 
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Nephanalysis  (3DNEPH)  program  was  developed  at  the  Air  Force  Global 
Weather  Central  (AFGWC)  to  incorporate  the  tremendous  quantity  of 
satellite  sensed  cloud  data  and  conventionally  sensed  meteorological 
parameters  into  a three-dimensional  cloud  model  of  the  atmosphere. 

Basic  to  the  design  of  the  3DNEPH  is  the  assumption  that  satellite 
information  is  available  for  its  data  base  in  a timely  manner.  However, 
in  the  event  that  satellite  data  is  not  available,  the  3DNEPH  has  the  cap- 
ability of  extrapolating  past  analysis  until  such  time  as  satellite  data 
does  become  available. 

The  3DNEPH  program  is  built  as  a series  of  input  processors. 
These  processors  include  the  surface  data  processor,  radiosonde  obser- 
vation (RAOB)  processor,  aircraft  data  processor,  manual  data  processor, 
decision  tree  processor,  satellite  video  data  processor,  satellite  in- 
frared data  processor,  final  processor,  forecast  processor,  verification 
processor,  and  display  processor.  Because  of  the  modular  nature  of  the 
3DNEPH  program,  processors  can  be  added  to  or  deleted  from  the  system 
with  a minimum  of  programming  problems.  Descriptions  and  functions  of 
each  of  the  processors  can  be  found  in  the  AFGWC  Technical  Memorandum 
71-2  by  Coburn  (1971). 

The  horizontal  resolution  of  the  3DNEPH  program  is  limited  by 
the  resolution,  and  mapping  and  gridding  accuracy  of  the  input  satellite 
data.  The  hemispheric  grid  chosen  for  the  3DNEPH,  which  was  compatible 
with  the  accuracy  of  its  input  satellite  data,  was  a 512  x 512  array  cen- 
tered at  the  north  (south)  pole  of  a polar  stereographic  map  and  having 
a distance  between  grid  points  of  40  km  at  60°  latitude.  This  512  x 512 


grid  was  further  subdivided  into  64  squares  (boxes)  so  that  the  finest 
mesh  was  4096  grid  points  over  each  hemisphere.  Once  again,  each  grid 
point  contains  information  representative  of  a 40  km  square  centered  at 
the  grid  point  when  at  60°  latitude.  The  vertical  resolution  of  the 
3DNEPH  program  divides  the  atmosphere  into  15  layers.  The  first  six  lay- 
ers are  terrain  following  layers  and  the  last  nine  layers  are  categorized 
in  feet  above  mean  sea  level  (MSL).  In  addition  to  the  information  given 
for  cloud  amounts  in  the  layers,  seven  additional  pieces  of  data  are  given 
at  each  point.  These  data  denote  information  about  the  cloud  types, 
maximum  tops  and  minimum  bases,  the  current  weather  and  the  total  cloud 
cover. 

The  use  of  the  3DNEPH  as  a source  of  comparison  to  the  technique 
developed  here  has  many  shortcomings.  These  problems  include  timeliness, 
layer  resolution,  dependence  on  subjective  observations,  and  the  horizon- 
tal spreading  of  the  data.  The  3DNEPH  was  used  as  a comparison,  however, 
since  no  other  known  source  containing  the  routine  analysis  of  cloud 
parameters  is  available. 

The  ever-changing  cloud  scene  at  a point  makes  the  timeliness  of 
the  3DNEPH  and  the  HIRS  data  very  important.  The  3DNEPH  is  produced  every 
three  hours,  and  if  there  is  no  new  data  for  an  update,  the  previous  an- 
alysis or  a cloud  forecast  is  used  for  the  new  analysis.  This  problem  is 
minimized  since  we  are  considering  overcast  cases  in  which  the  cloud  scene 
at  a point  will  probably  not  change  significantly  over  a short  period  of 
time. 

A second  problem  area  is  the  layer  resolution  of  the  3DNEPH. 
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The  resolution  can  be  as  poor  as  4000  feet  for  middle  cloud  layers  and 
as  thick  as  20,000  feet  for  high  cloud  layers.  If  a thin  layer  of  cloud 
is  reported  to  be  in  one  of  the  layers,  the  model  must  automatically 
fill  the  entire  layer  with  clouds.  As  noted  earlier,  an  error  of  4,000 
ft  (1.2  km)  is  one-half  of  the  thickness  of  middle  cloud  that  can  be  de- 
tected by  the  theoretical  model.  In  the  case  of  cirrus  clouds,  a 20,000 
ft  error  is  equivalent  to  6.1  km,  a greater  thickness  than  the  model  can 
handle.  This  large  error  can  be  somewhat  offset  when  the  maximum  cloud 
top  parameter  is  considered.  At  middle  and  high  cloud  top  levels,  the 
coding  for  this  parameter  is  accurate  to  5,000  ft  or  1.5  km.  For  lower 
middle  cloud  tops,  this  parameter  is  accurate  to  1,000  ft  or  0.3  km. 

The  use  of  subjective  weather  observations  plays  an  important 
role  in  the  accuracy  of  the  3DNEPH  clouds.  If  there  is  no  current 
satellite  data,  the  cloud  top  cannot  be  determined  and  persistence 
prevails,  and  the  cloud  thickness  is  over  or  under-estimated.  Another 
case  would  be  when  satellite  data  is  available  but  no  surface  observation 
is  present  (over  oceans).  Here,  the  satellite  can  determine  the  cloud  top 
height  only,  and  no  information  about  the  cloud  depth  can  be  determined. 

Another  problem  is  that  surface  observations  are  spread  horizon- 
tally. The  3DNEPH  grid  point  represents  a 40  km  box  at  60°N.  If  there 
is  no  current  satellite  data  and  the  observation  network  is  not  at  a 40 
km  resolution,  some  influence  of  the  surface  observation  is  spread  to 
adjacent  points  in  the  3DNEPH.  This  can  be  seen  by  examining  an  area 
of  3DNEPH  data  where  the  cloud  thickness  will  be  constant  over  a large 
area.  Steps  taken  to  minimize  these  problem  areas  will  be  discussed 
below. 
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4.4.2  High  resolution  infrared  sounder  data.  HIRS  data  is  routinely 
processed  at  NOAA/NESS  and  is  available  on  nine-track,  1600  bpi  tapes. 

The  data  is  packed  in  such  a way  that  six  or  seven  orbits  are  available 
on  each  tape.  These  tapes  contain  located  and  calibrated  radiance  values 
for  all  17  channels.  The  Nimbus  VI  User's  Guide  gives  a detailed  explan- 
ation of  the  calibration  procedure  used  on  the  radiance  measurements  in 
each  channel.  The  resolution  of  the  HIRS  data  decreases  with  increasing 
scan  angle.  This  decrease  in  resolution  of  a scan  spot  is  from  approxi- 
mately 23  km  at  nadir  to  31  km  at  a scan  angle  of  36.9°. 

The  data  set  used  in  this  analysis  was  for  20-30  August  1975, 
a time  when  all  channels  of  the  HIRS  instrument  were  operating  properly. 
Further,  it  was  required  to  have  data  over  North  America  since  the  com- 
parison was  being  used  with  the  3DNEPH  which  depends  heavily  upon  good 
surface  observations.  We  therefore  requested  data  from  20°-55°  N and 
80°-150°  W which  includes  both  land  and  ocean  areas.  To  further  com- 
plicate the  data  selection,  examination  of  satellite  photographs  indi- 
cated a general  lack  of  middle  and  high  clouds  over  this  region  for  this 
time  period  except  in  areas  north  of  40°  N.  As  a result,  all  of  the  data 
used  in  the  analysis  are  along  a storm  track  north  of  40°  N. 

A total  of  26  passes  were  provided  in  the  data  set  of  which 
five  were  analyzed  for  cloud  information  at  the  two  scan  angles  closest 
to  nadir.  To  check  consistency  of  the  data  from  channel  to  channel 
within  a pass,  the  five  passes  used  in  the  analysis  were  displayed  by 
using  an  overprinting  technique  to  simulate  grey  shading  on  the  line 
printer.  All  the  passes  used  in  subsequent  analyses  were  overprinted 
to  verify  the  HIRS  data  by  comparing  the  overprints  with  surface  analysis 
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The  basic  criteria  for  the  selection  of  3DNEPH  data  for  compari- 
son with  the  HIRS  data  was  the  reliability  of  the  3DNEPH  data.  This 
criteria  dictated  that  data  be  chosen  where  the  3DNEPH  had  good  surface 
data  as  well  as  satellite  data.  Additionally,  the  areas  of  interest 
must  be  in  the  presence  of  synoptic  scale  weather  so  that  continuous 
cloud  decks  would  be  insured.  The  use  of  HIRS  data  points  at  or  near 
nadir  would  insure  the  least  amount  of  error  in  applying  the  real  data 
to  the  theoretical  calculations,  since  the  fitting  of  the  theoretical 
upwelling  radiances  to  different  scan  angles  would  be  most  accurate  at 
nadir. 

Based  on  these  criteria,  five  passes  on  five  consecutive  days 
were  chosen  for  this  comparison.  These  days  included  August  21  through 
August  25,  1975.  The  pass  on  each  day  that  included  Western  North  America 
and  was  closest  to  a 3DNEPH  analysis  time  was  used.  Nimbus  VI  has  an 
approximate  equator  crossing  at  1000  local  time  on  its  ascending  pass  and 
corresponds  to  an  1800Z  (GMT)  analysis  time  for  the  3DNEPH.  Box  number 
44  of  the  3DNEPH  covers  this  area  and  has  the  following  latitudes  and 
longitudes  for  the  upper  left,  upper  right,  lower  left  and  lower  right, 
respectively,  52°N,  125°W,  62°N,  80°W,  31 °N,  106°W,  and  36°N,  80°W. 
Inspection  of  satellite  pictures  on  these  five  days  show  that  all  the 
clouds  are  north  of  40°N  and  a majority  of  them  are  north  of  45°N. 

Hence,  all  the  points  used  in  the  comparison  are  at  these  latitudes. 

A clear  column  radiance  representative  of  the  pertinent  area 
was  chosen  using  the  satellite  pictures  in  conjunction  with  the  HIRS 
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data  for  each  day  under  investigation.  These  clear  column  radiances  were 
used  to  reduce  the  cloudy  radiances  to  ratios  as  was  discussed  in  a pre- 
vious section.  Note  that  both  theoretical  calculations  and  the  real  data 
are  for  mid-latitude  summer  conditions. 

The  colocation  of  the  HIRS  data  with  the  3DNEPH  was  accomplished 
using  the  latitude  and  longitude  from  the  HIRS  data.  These  quantities 
were  then  converted  to  3DNEPH  grid  location.  Since  this  location  is  never 
an  exact  3DNEPH  point,  the  3DNEPH  point  used  for  the  comparison  was  that 
point  closest  to  the  calculated  grid  location. 


Cloud  Type  and  Mass  Determination  from  HIRS  Data 


4.5.1  Cloud  type.  All  of  the  HIRS  data  applied  to  the  theoretically 
based  empirical  relationships  are  from  the  cloudy  areas,  north  of  40°N. 

The  slope  and  y-intercept  of  each  case  were  calculated.  The  cases  with  the 
magnitude  of  the  slope  greater  than  the  magnitude  of  the  y-intercept  were 
assumed  to  be  cirrus  clouds.  The  cases  with  the  magnitude  of  the  y-inter- 
cept greater  than  the  magnitude  of  the  slope  were  assumed  to  be  middle  or 
lower  clouds.  These  cases  were  then  examined  in  conjunction  with  the  IR 
pictures.  In  all  cases,  the  empirically  derived  relationships  that  were 
based  on  the  theoretical  model  indicated  the  proper  cloud  type. 

The  theoretical  empirical  model  indicated  the  presence  of  middle 
clouds  in  the  vicinity  of  52°N,  115°W  to  43°N,  100°W  on  August  21.  Exam- 
ination of  the  IR  picture  not  shown  here  verifies  the  presence  of  lower 
clouds  in  this  area.  For  August  22,  the  model  indicated  a mixture  of  cirrus 
and  lower  clouds  from  52°N,  104°W  to  49°N,  102°W  and  then  extensive  lower 
clouds  to  43°N,  99°W.  This  same  pattern  can  be  noted  by  examining  Figure  4.7. 
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In  the  vicinity  of  50°N  the  IR  picture  shows  an  area  of  some  cirrus  clouds 
and  some  warmer  lower  clouds.  South  of  49°N  there  are  only  lower  clouds  in 
the  IR  picture  of  Figure  4.7.  It  also  indicates  little  cloud  south  of  43°N. 
One  of  the  clear  points  in  this  area  was  used  for  the  clear  column  radiance 
in  the  ratioing  for  this  day.  The  IR  picture  for  August  23  (not  shown  here) 
indicates  extensive  cirrus  on  the  subtrack  of  the  satellite  from  53°N,  120°W 
to  42°N,  115°W.  The  model  shows  all  cirrus  clouds  in  this  area  except  for 
the  indication  of  some  middle  clouds  near  53°N.  Further  examination  of  the 
IR  picture  on  this  day  reveals  some  breaks  in  the  cirrus  cloud  in  this 
area.  On  August  24,  cirrus  was  again  indicated  from  52°N,  109°W  to  48°W, 

1 06°W  and  middle  clouds  were  present  south  of  47°N.  Examination  of  the 
satellite  picture  for  that  day  not  shown  here  verifies  this  analysis.  The 
clear  column  radiance  for  this  day  was  taken  from  an  area  just  south  of 
47°N  along  the  subtrack  of  the  satellite.  The  model  analyzed  middle  clouds 
from  53°N,  97.5°W  to  47°N,  94.5°W  and  cirrus  clouds  from  44°N,  93°W  to  3d°N, 
91.7°W  for  August  25.  This  again  is  represented  well  on  the  IR  picture  of 
Figure  4.3.  The  clear  column  radiance  for  this  day  was  from  near  38°N 
and  9I°W.  The  ice  and  liquid  water  content  for  this  pass  along  with  the 
pass  for  August  22  are  mapped  in  a subsequent  section. 

4.2  Cloud  mass.  To  determine  the  cloud  mass  the  data  described  in 

the  previous  section  was  divided  into  middle  and  high  clouds.  The  3DNFPH 
points  that  correspond  to  these  data  cases  were  then  selected  from  Box 
44  of  the  3DNEPH  at  1800Z  for  the  five  days  analyzed.  In  all  cases  the 
3DNEPH  indicated  lower  clouds  below  the  cirrus  or  middle  cloud.  This 
indicates  one  of  the  basic  problems  discussed  earlier  about  the  3DNFPH, 
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that  with  only  surface  analysis  the  cloud  top  is  over  estimated  and  with 
only  satellite  data  the  cloud  base  is  lower  than  it  should  be,  so  that 
the  error  can  be  greater  than  the  thickness  of  the  cloud.  It  is  virtually 
impossible  to  correct  this  problem  without  instrumentation  that  can  see 
through  clouds.  In  addition,  the  spreading  of  an  observation  to  many 
3DNEPH  points  tends  to  smooth  the  area  out  so  that  cloud  thicknesses  are 
constant  over  large  areas.  To  minimize  these  problems,  the  30NEPH  cloud 
thickness  data  were  used  as  either  cirrus  or  middle  cloud.  As  noted 
earlier,  this  would  not  affect  the  ice  and  water  content  calculation 
significantly  because  middle  cloud  in  the  presence  of  thin  cirrus  has 
the  same  appearance  as  moderately  thick  cirrus.  The  effect  of  lower 
clouds  on  middle  cloud  ratios  remained  virtually  unchanged  when  comparing 
the  slopes  and  y-intercepts  derived  from  the  empirical -theoretical  approach. 

The  3DNEPH  thicknesses  were  parameterized  by  assuming  that  the 

actual  base  of  the  cloud  corresponded  to  the  theoretical  base  of  the  cloud. 

This  was  3.1  km  for  middle  cloud  and  6.7  km  for  cirrus  clouds.  The  cloud 

thickness  in  the  middle  cloud  case  was  then  obtained  by  adding  the  terrain 

height  to  the  middle  cloud  base  and  subtracting  this  from  the  layer  top 

or  the  maximum  cloud  top  parameter.  This  was  then  multiplied  by  1.5  x 
2 -2  -1 

10  gm  m km  to  obtain  the  liquid  water  content.  The  cirrus  cloud  thick- 
ness was  obtained  by  subtracting  6.7  km  from  the  height  of  the  top  of  the 
layer  containing  clouds  or  the  maximum  top  parameter  whichever  was  smaller. 

O -1 

This  result  was  then  multiplied  by  28.3  gm  m km  to  obtain  the  ice 
content.  This  should  minimize  as  much  as  possible  the  thickness  grouping 
in  the  3DNEPH. 

Figure  4.8  is  a graph  comparing  the  model  ice  and  liquid  water 
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contents  with  the  3DNEPH  derived  ice  and  water  contents.  Part  a of  the 
graph  is  for  middle  cloud  cases  and  part  h is  for  cirrus  cloud  cases. 

In  both  parts,  both  the  3DNEPH  and  the  theoretical-empirical  relation- 
ships indicated  middle  or  cirrus  clouds.  The  ordinate  is  for  the  3DNEPH 
and  is  labeled  in  both  thickness  and  equivalent  ice  or  liquid  water  con- 
tent based  on  the  ice  and  liquid  water  contents  used  in  the  theoretical 
calculations.  The  abscissa  is  for  the  empirical -theoretical  calculations 
and  is  also  labeled  in  thickness  and  equiyalent  liquid  water  content  for 
part  a and  equivalent  ice  content  for  part  b. 

Comparison  of  the  3DNEPH  derived  water  contents  in  Figure  4.R 
with  the  empirically  derived  theoretical  water  contents  is  somewhat  vari- 
able. Although  the  linear  relationship  desired  is  not  present  some 
correlation  can  be  noted.  The  data  points  in  the  comparison  primarily 
occur  at  four  3DNFPH  thicknessess  while  the  model  derived  water  contents 
covered  a wide  range.  This  grouping  is  present  since  the  3DNEPH  has 
only  four  middle  cloud  layers  in  its  data  base.  The  further  parameteri za- 
tions  done  to  the  3DNEPH  in  this  analysis  mads  ;t  possible  to  increase 
the  number  of  water  contents  that  could  be  obtained  from  the  3DNEPH . The 

error  range  of  the  comparisons  for  part  a of  Figure  4.8  is  up  to  2.25  x 
2 -2 

10  gm  m . Analysis  of  Figure  4.8  also  indicates  the  grouping  of  the 
30NEPH  data  at  basically  three  ice  content  values  although  the  range 
from  the  empirical-theoretical  calculations  vary  widely  over  a given 
30NEPM  value.  The  error  range  in  Figure  4.8  is  as  much  as  84.8  gm  m . 

We  note  that  there  appears  to  be  some  agreement  between  the  two  independent 
methods  of  computing  cloud  ice  or  liquid  water  content. 

The  above  comparisons  point  out  the  simple  fact  that  new  methods 
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Figure  4.8  Comparison  of  3DNEPH  middle  cloud  water  content 
with  empirical-theoretical  middle  cloud  water 
content  from  HIRS  in  Part  a and  comparison  of 
3DNEPH  cirrus  cloud  ice  content  with  empirical- 
theoretical  cirrus  ice  content  from  HIRS  in 
Part  b. 


for  determining  cloud  structures  from  satellite  data  cannot  be  properly 
verified  with  the  current  data  bases  that  are  available.  The  cloud  tvninn 
derived  in  this  analysis  can  be  compared  with  simultaneous  cloud  photo- 
graphs and  these  comparisons  appear  to  be  good  in  light  of  discussions 
in  the  previous  section.  The  verification  of  ice  and  liquid  water  content 
will  require  well  planned  and  controlled  measurements  of  cloud  structures 
within  the  satellite  field-of-view.  With  this  type  of  verification,  it 
will  be  possible  to  reformulate  parameterizations  of  theoretical  calcula- 
tions to  improve  the  accuracy  of  these  estimates  of  cloud  compositions. 

Until  this  can  be  accomplished,  the  use  of  theoretical -empirical  relation- 
ships should  prove  to  be  more  accurate  than  purely  statistical  methods. 

4.5.3  Cloud  moisture  mapping.  The  use  of  the  parameterizations  described 
in  the  previous  section  and  applied  to  HIRS  data  in  this  section  can  be  used 
to  map  ice  and  water  content  of  middle  and  cirrus  clouds.  To  perform  this 
mapping  on  a global  scale  over  long  periods  of  time  would  be  the  desired 
result.  However,  to  accomplish  this  on  a routine  basis  would  reauire  the 
availability  of  global  satellite  data  and  a large  block  of  dedicated  com- 
puter time.  In  this  section,  portions  of  two  orbits  over  central  North 
America  are  mapped  to  illustrate  the  concept. 

The  maps  used  for  these  projects  are  polar  stereographic  and 
it  is  therefore  possible  to  compare  rather  easily  to  NDAA  4 satellite 
pictures.  In  each  of  the  cases  the  clear  column  radiance  at  satellite 
nadir  was  used  to  perform  the  ratioing  over  ten  scan  spots  on  each  side 
of  nadir.  Beyond  these  scan  angles  a clear  column  radiance  represen- 
tative of  a large  scan  angle  should  be  used  to  perform  the  ratioinq  to 
derive  the  amount  of  liquid  water  or  ice  content. 
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As  noted  earlier,  the  slope  and  y-intercept  of  the  linear  fit 
determine  whether  middle  or  high  clouds  are  present.  If  the  magnitude 
of  the  slope  is  greater  than  the  y-intercept  cirrus  cloud  is  assumed. 

A further  criterion  established  for  mapping  was  to  eliminate  data  points 
where  the  magnitude  of  the  slope  was  greater  than  25,  or  the  correlation 
coefficient  was  less  than  0.70.  After  comparing  many  cases  of  actual  data 
with  theoretical  calculations,  we  found  that  cases  that  do  not  meet  the 
above  criterion  were  low  cloud  which  were  not  addressed  in  this  study. 

For  passes  on  August  22  and  August  25,  the  latitudes  and  lonai- 
tudes  of  the  data  points  that  met  the  above  criterion  were  plotted  on  the 
polar  projection  and  then  analyzed.  In  both  cases  the  cirrus  and  middle 
cloud  are  mapped  on  the  same  figure.  The  isolines  are  labeled  with  verti- 
cal ice  or  liquid  water  content  based  on  the  constants  derived  from  para- 
meterizations  of  the  theoretical  calculations.  The  solid  lines  are  isolines 
of  middle  cloud,  dashed  lines  are  for  cirrus  cloud  and  the  darker  lines  are 
separation  of  middle  and  high  cloud  analysis. 

The  analysis  for  August  22  shows  two  areas  of  overcast  conditions 
under  the  satellite.  The  northern  area  is  over  southern  Canada  and  extreme 
northern  Montana  and  North  Dakota  and  is  characterized  by  both  middle  and 
high  clouds.  Since  the  data  is  cut  off  to  the  north,  east  and  west,  the 
analysis  shows  sharp  break  points  in  these  directions.  NDAA  4 pictures 
indicate  that  the  southern  cloud  areas  analyzed  for  this  pass  have  little 
cirrus  cloud  present  which  is  in  agreement  with  the  present  analysis.  The 
theoretical  analysis  also  shows  middle  cloud  less  than  1 km  thick  both  south 
of  the  northern  area  and  north  of  the  southern  area  of  clouds. 

Analysis  of  the  figure  for  August  25  shows  a similar  set  of  two 
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Figure  4.9  Cloud  ice  and  water  content  based  on  empirical- 
theoretical  calculations  for  August  22,  1975. 
Dotted  lines  are  for  cirrus  clouds  (gm  m*2)  and 
solid  lines  are  for  middle  clouds  (gm  m-2).  The 
heavy  line  is  the  cloud  boundaries. 


Figure  4.10  Cloud  ice  and  water  content  based  on  empirical  - 
theoretical  calculations  for  August  25,  1975. 
Dotted  lines  are  for  cirrus  clouds  (gm  m"z)  and 
solid  lines  are  for  middle  clouds  (gm  m*2).  The 
heavy  line  is  the  cloud  boundaries. 
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distinct  cloud  areas.  The  southern  area  is  the  cirrus  on  top  of  an  active 
area  of  thunderstorms  and  the  sharp  southern  edge  of  this  cloud  mass  cor- 
responds very  well  to  the  NOAA  4 picture.  The  northern  cloud  mass  is 
dominated  by  middle  cloud  with  some  cirrus  on  the  southwest  side.  This 
again  agrees  well  with  the  NOAA  4 satellite  picture.  On  the  southern 
edge  of  the  northern  cloud  mass  the  middle  cloud  thicknessess  fall  off 
very  rapidly  on  the  picture  and  the  mapping. 

With  proper  verification,  the  use  of  this  technique  on  a global 
scale  is  feasible.  Proper  verification,  as  noted  in  the  previous  section, 
would  require  the  use  of  cloud  physics  measurements  in  the  satellite  field- 
of-view.  These  types  of  measurements  would  be  required  in  a variety  of 
atmospheric  conditions.  The  method  presented  here  seems  to  lend  itself  to 
real  time  analysis  since  the  computer  time  requirement  has  been  minimized 
by  the  empirical  parameterizations.  In  addition  the  ratioing  has  decreased 
the  effect  of  the  real  atmosphere  on  the  cloud  moisture  determination.  The 
effects  of  middle  clouds  with  colder  temperatures  than  those  used  in  the 
theoretical  calculations  are  taken  into  account  since  colder  clouds  repre- 
sent clouds  thicker  than  the  model  can  handle  and  middle  clouds  usually 
build  from  the  base.  It  should  be  noted  that  only  the  climatological  mid- 
latitude summer  atmospheric  profiles  were  used  in  theoretical  analyses. 

The  technique  developed  here  when  compared  to  microwave  tech- 
niques has  both  advantages  and  disadvantages.  Two  advantages  are  the 
higher  resolution  of  the  HIRS  instrument  when  compared  to  the  Nimbus 
VI  SCAMS  instrument  and  the  detection  of  cirrus  clouds  that  are  trans- 
parent in  the  microwave  portion  of  the  spectrum.  Another  advantage  of 
this  technique  is  that  it  can  be  applied  over  both  land  and  oceans,  while 
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microwave  interpretation  of  clouds  is  restricted  to  oceans  because  of 
the  change  of  land  emissivity.  Disadvantages  of  this  technique  when 
compared  to  microwave  include  the  opacity  of  water  clouds  in  the  infra- 
red portion  of  the  spectrum  and  the  sensitivity  of  the  technique  to  the 
water  cloud  temperature.  By  combining  infrared  and  microwave  measurements 
along  with  the  proper  cloud  physics  verification  data,  operationally 
significant  recovery  of  cloud  parameters  from  satellite  data  appears 
feasible. 

4.6  Conclusions 

A theoretical  model  that  calculates  transfer  of  spectral  infra- 
red radiation  based  on  the  discrete-ordinate  method  developed  in  Section 
2 was  modified  to  include  cirrus  and  middle  clouds  and  absorbing  gases. 
This  analysis  used  one  model  atmosphere  and  one  middle  and  high  cloud 
type  for  all  calculations.  The  top  of  the  cirrus  cloud  and  the  base  of 
the  middle  cloud  were  held  constant  to  minimize  computer  time.  The 
layers  in  the  transfer  computations  were  assumed  to  have  constant  temper- 
atures so  that  the  effects  of  the  change  in  cloud  thickness  on  the  up- 
welling  radiance  could  be  studied.  Various  thicknesses  and  combinations 
of  middle  and  cirrus  clouds  were  used  in  the  analysis.  This  model  was 
then  applied  to  the  HIRS  channels  of  the  Nimbus  VI  satellite.  A case 
study  using  the  actual  HIRS  radiance  data  shows  that  it  is  feasible  to 
estimate  the  cirrus  cloud  thickness  employing  the  combination  of  infrared 
channels  with  the  assistance  of  IR  and  visible  cloud  pictures  from  the 
NOAA  IV  satellite. 

The  theoretical  calculations  of  upwelling  radiance  were  divid- 
ed by  their  clear  column  radiances  to  give  a ratio  for  each  channel  in 
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which  the  atmospheric  temperature  effect  was  minimized.  These  ratios 
were  then  fitted  with  a straight  line  for  each  theoretical  cloud  thick- 
ness for  both  cirrus  and  middle  clouds.  It  was  found  that  cases  where 
the  slope  was  greater  than  the  magnitude  of  the  y- intercept  derived  from 
transfer  calculations  corresponded  to  the  cirrus  clouds.  The  reverse 
of  the  slope  and  y-intercept  relation  corresponded  to  middle  clouds  or 
lower  clouds.  The  resulting  slopes  and  y-intercepts  for  different  cloud 
thicknesses  for  the  two  cloud  types  were  then  fitted  with  a logarithmic 
function.  This  made  it  possible  to  infer  ice  and  water  content  in  the 
case  of  high  and  middle  cloud*:,  respectively,  for  a given  thickness  based 
on  the  theoretical  calculations. 

These  theoretical  parameterizations  were  applied  to  five  days  of 
HIRS  data.  The  empirical  parameterization  method  appears  to  be  successful 
in  identifying  high  and  middle  clouds.  The  resulting  cloud  ice  and  liquid 
water  contents  were  compared  to  the  3DNEPH  to  indicate  whether  this  tech- 
nique could  be  applied  in  an  operational  mode.  It  was  found  that  the 
comparisons  with  the  3DNEPH  were  marginal.  This  can  be  traced  to  the 
3DNEPH  and  the  methods  used  to  determine  cloud  thicknesses  as  well  as 
approximations  in  the  theoretical  calculations  and  the  parameterization 
of  these  calculations.  Two  days  of  ice  and  water  content  for  actual 
HIRS  data  were  then  mapped.  These  results  were  good  when  comparing  the 
resulting  maps  with  the  corresponding  NOAA  IV  infrared  and  visible  satel- 
lite pictures. 

The  technique  developed  in  this  report  is  ideally  suited  to 
enhance  the  cloud  parameterizations  in  the  3DNEPH.  As  noted  earlier 
the  3DNEPH  is  a modular  program  that  processes  all  kinds  of  meteorological 
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data.  The  only  satellite  data  that  is  currently  used  is  broad  band  visible 
and  infrared  channels  of  a NOAA  4 type  scanning  radiometer.  The  addition 
of  this  technique  would  use  satellite  data  that  is  not  in  current  use 
for  cloud  detection  and  it  would  give  valuable  input  into  the  final  re- 
sults. Since  the  3DNEPH  is  a routinely  produced  product,  much  of  the 
software  required  to  implement  this  procedure  is  already  in  existence 
and  the  operational  development  could  be  carried  out  independently  of 
3DNEPH  production.  As  noted  earlier  the  resolutions  of  the  HIRS  and  the 
3DNEPH  are  similar  so  that  in  an  operational  implementation  the  3DNEPH 
could  be  used  to  further  verify  the  technique. 

Although  an  objective  way  of  deriving  cloud  type  information  and 
cloud  ice  or  water  content  has  been  illustrated  in  this  analysis,  a more 
reliable  method  of  verification  must  be  found.  However,  verification  of 
satellite  sensing  techniques  requires  carefully  designed  field  experiments 
in  which  highly  reliable  cloud  parameters  could  be  obtained  under  the 
satellite  pass.  With  sufficient  cases  from  which  cloud  parameters  may 
be  derived  locally,  intercomparisons  with  satellite  derived  values  may 
be  carried  out  to  establish  the  statistical  significance  of  this  techni- 
que to  parameterize  satellite  measurements. 
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ABSTRACT 

A retrieval  technique  is  presented  for  the  determination  of  the  surface  temperature,  the  thickness  and 
transmissivity  of  cirrus  clouds,  and  the  fraction  of  the  cirrus  cloudiness  by  means  of  four  observed  upwelling 
radiances  in  the  10  pm  window  region.  On  the  basis  of  radiative  transfer  calculations  for  mean  wavenumbers 
of  900,  950,  1 100  and  1 150  cm-1,  assumptions  are  made  in  the  theoretical  retrieval  analyses  that  water  vapor 
effects  al>ove  cirrus  clouds  are  negligible  and  that  ratios  of  the  transmissivities  are  linear  functions  of  the 
cloud  thickness.  Krror  analyses  employing  climatological  data  reveal  that  independent  random  errors  in 
temperature  and  humidity  profiles  introduce  insignificant  errors  in  the  four  resulting  parameters.  The  result- 
ing errors  caused  by  random  errors  in  the  expected  upwelling  radiances,  however,  depend  upon  their  standard 
deviations.  Once  the  thickness  and  the  transmissivity  at  a given  wavenumber  of  a cirrus  cloud  have  been 
determined,  we  illustrate  that  the  vertical  ice  content  may  be  estimated  assuming  that  ice  particles  are 
randomly  oriented  in  a horizontal  plane. 


1.  Introduction 

Cirrus  clouds  have  been  noted  to  introduce  serious 
difficulties  in  remote  sensing  of  atmospheric  tempera- 
ture and  humidity  profiles  and  surface  conditions, 
owing  to  their  semi-transparent  appearance  in  the 
visible  as  well  as  in  the  infrared  region.  It  is  extremely 
important,  therefore,  to  derive  the  reliable  and  ac- 
curate radiation  properties  of  cirrus  clouds  from  an 
independent  set  of  sounding  frequencies  in  conjunc- 
tion with  the  atmospheric  parameter  evaluation. 
Moreover,  determination  of  the  vertical  ice  content 
over  the  global  atmosphere  is  equally  important  from 
the  point  of  view  of  climatology  studies  and,  perhaps, 
numerical  weather  prediction. 

The  paper  by  Houghton  and  Hunt  (1971)  ap- 
parently was  the  first  one  to  explore  the  passive 
remote  sensing  possibility  of  ice  clouds  by  means  of 
two  wavelengths  in  the  far  infrared.  Liou  (1974) 
discussed  emission  and  transmission  properties  of  cirrus 
clouds  in  the  10  nm  window  region  in  conjunction 
with  the  remote  sensing  potential  from  satellites. 
Bunting  and  Conover  (1974)  proposed  a simple  means 
for  the  estimation  of  the  vertical  ice  content  of  cirrus 
clouds  assuming  the  exponential  attenuation  of  IR 
radiation. 

With  respect  to  the  atmospheric  parameter  deter- 
mination, Chahine  (1974)  presented  a numerical 
method  to  derive  vertical  temperature  profiles  in 
cloudy  atmospheres  by  means  of  radiances  obtained 
from  two  partially  overlapping  fields  of  view.  Taylor 
(1974)  described  an  approach  employing  soundings  at 


two  different  zenith  angles  for  the  determination  of 
temperature  profiles  in  the  presence  of  clouds.  Both 
of  these  two  important  studies  concerning  the  remote 
sensing  of  cloudy  atmospheres  involved  attempts  to 
remove  the  cloud  effect  rather  than  to  incorporate 
it  into  the  radiative  transfer  equation. 

Cirrus  clouds  consist  of  non-spherical  ice  crystals 
of  various  sizes,  possibly  randomly  oriented  in  a hori- 
zontal plane.  Other  unknown  variables  include  the 
ice-crystal  concentration,  the  cloud  thickness  and  the 
location  of  the  cloud  in  the  atmosphere.  Moreover, 
for  a given  upwelling  radiance  at  the  satellite  point 
of  view',  there  will  be  additional  unknown  variables 
associated  with  the  atmospheric  temperature  and 
gaseous  profiles.  In  view  of  a large  number  of  un- 
known variables,  parameterization  equations  describing 
the  radiation  field  of  an  atmosphere  containing  cirrus 
clouds  have  to  be  formulated.  Retrieval  of  the  un- 
known parameters,  each  of  which  may  represent  a 
combination  of  several  variables,  may  be  carried  out 
from  satellite  radiance  observations. 

In  this  paper,  I would  like  to  present  some  hypo- 
thetical analyses  and  calculations  concerning  the  de- 
termination of  the  structure  and  composition  of  cirrus 
clouds  from  a set  of  synthetic  radiance  observations 
in  the  window  region. 

2.  Theoretical  analyses 

We  assume  that  within  tne  field  of  view'  of  the 
satellite  radiometer  the  atmosphere  contains  >)  portion 
of  cloudiness.  The  monochromatic  upwelling  radiance 
measured  by  the  satellite  radiometer  at  the  top  of  a 
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Fic.  1.  Values  of  the  integral  term  and  transmission  function 
shown  in  Eq.  (2).  The  integral  term  is  evaluated  from  the 
surface  to  the  level  s,  while  the  transmission  function  of  water 
vapor  is  calculated  from  the  top  of  the  atmosphere  to  the  level  t. 
The  calculations  are  carried  out  for  wavenumbers  of  850,  950, 
1150  cm"1  employing  a tropical  atmosphere. 

partially  cloudy  atmosphere  is  given  by 

7r(«)=n/j(«)+(i-fl)/r(»),  (i) 

where  the  upwelling  radiance  arising  from  the  clear 
portion  of  the  atmosphere  is  simply 

/?(»)  = /U7'.M»,0)+  \ /L[r(2)}/T,Ks).  (2) 

J 1-0 

And  the  upwelling  radiance  due  to  the  cloudy  portion 
of  the  atmosphere  may  be  expressed  as 

/;(»)  = r.7,(sk)r,(=c,t,)+  f (3) 

J —t, 

In  this  equation  the  cloud  transmissivity  t[  is  defined 
as  the  ratio  of  the  upwelling  radiance  at  the  cloud 
top,  i.e.,  I. (it),  to  that  at  the  cloud  base,  and  it  can 
be  written  as 

r‘.  = /.(*,)//,(**),  (4) 

where  and  z,  denote  the  cloud-base  and  cloud-top 
heights,  respectively,  and  the  upwelling  radiance 
reaching  the  cloud  base  is  given  by 

M*»)-W)r.(*»,0)+  [ ‘ PlT(t »,(*.,*),  (5) 

J —t 


In  Eqs.  (3)-( 5),  T,  is  the  surface  temperature, 
t,  the  transmission  function  of  gases  which  will  be 
discussed  later,  and  the  Planck  function  in  the 
wavenumber  domain  is 

B,(T)  = aS/ (e'’”- 1),  (6) 

with  a=  1.1272X10“*  erg  cm*  s-1  and  fr=  1.4389  cm  K. 
It  should  be  noted  that  /,(z()  in  Eqs.  (3)  and  (4) 
must  be  evaluated  from  the  transfer  equation  including 
absorption,  emission  and  scattering  of  gases  and  cloud 
particles  subject  to  the  radiation  boundary  conditions 
of  the  top  and  base  of  the  cloud  concerned  (Liou, 
1974). 

Cirrus  clouds  are  normally  fairly  high  in  the  atmo- 
sphere with  their  top  heights  on  the  order  of  It)  km 
or  higher.  Thus,  we  may  select  spectral  channels  in 
the  window  region  where  the  effect  of  water  vapor 
absorption  above  the  cirrus  cloud  layer  can  be  ne- 
glected such  that 

/•*“*> 

j I — tf  "*  (7) 

T,(«,z,)«l. 

Justification  of  Eq.  (7)  is  illustrated  in  Fig.  1 (sec 
also  Section  3a  for  water  vapor  absorption).  We  have 
employed  a tropical  atmospheric  profile  (McClalchey 
et  al.,  1971)  in  the  calculations.  The  upper  diagram 
represents  the  results  of  the  first  equation  in  (7)  for 
the  integration  limits  from  the  surface  to  a level  z. 
Three  wavenumbers  of  850,  950  and  1150  cm-1  were 
used.  It  is  evident  that  the  integral  term  in  Kq.  (2) 
has  negligible  contribution  above  a height  of  about 
5 km,  which  is  normally  lower  than  the  base  heights 
of  cirrus  clouds.  The  lower  diagram  represents  the 
transmission  term  in  Eq.  (2).  Clearly,  we  see  that 
above  about  8 km,  the  transmission  function  t,  of 
water  vapor  for  a moist  atmosphere  is  very  close  to  1. 
Consequently,  we  would  expect  that  our  first  assump- 
tions can  be  applied  to  other  atmospheric  profiles 
without  difficulties. 

Thus  the  observed  upwelling  radiance  at  the  top 
of  a partially  cloudy  atmosphere  denoted  in  Eq.  (1) 
becomes 

7r(«)  ~ (1  - n+  vri)  ( B.(T.) T,(ih,(I) 

+ \ HlT(t)y  r.(i^)\,  (8) 

J »-P 

where  the  unknown  parameters  are  the  fraction  of 
cloudiness,  the  surface  temperature,  the  cloud  trans- 
missivity which  is  wavenumber  dependent,  and  the 
tem|ierature  profile  and  the  transmission  function  of 
water  vapor.  Eq.  (8)  expresses  the  upwelling  radiance 
for  a single  wavenumber.  However,  a satellite  instru- 
ment can  distinguish  only  finite  bandwidths  At'(t>i,i2) 
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with  the  instrumental  slit  function  0(ri,Rj),  which  we 
shall  ignore  in  the  following  theoretical  development. 
The  finite  bandwidth  is  normally  so  small  that  the 
variation  of  Planck  function  with  respect  to  the 
wavenumber  can  be  neglected.  This  is  particularly 
evident  in  the  10  pm  window  region.  Thus,  the  up- 
welling  spectral  radiance  may  be  written  as 

/£(*>)  = [ 7r  (<■■>)—  = (1  - 9+  1^,) ( lh.{T.)rA,(zbt0) 

J ..  Ar 


+ 


(9) 


Here  we  also  note  that  the  variation  of  cloud  trans- 
missivity within  a small  spectral  interval  is  relatively 
small  in  view  of  the  slow  varying  refractive  indices 
of  ice  in  the  window  region.  Consequently,  it  is  physi- 
cally reasonable  to  use  a mean  wavenumber  to  cal- 
culate the  spectral  cloud  transmissivity. 

If  prior  knowledge  of  the  temperature  and  humidity 
profiles  in  the  atmosphere  were  available,  say  from 
climatological  data,  then  r,j  ,(**,())  and  the  integral 
term  in  Eq.  (9)  could  be  evaluated.  With  a measure- 
ment of  the  upwelling  radiance  from  the  satellite 
radiometer  in  the  window  channel,  there  are  only 
three  unknown  parameters,  tj,  t^,  and  T,.  How  many 
radiance  measurements  are  required  to  determine  all 
these  variables? 

We  shall  now  examine  the  cloud  transmissivity 
defined  in  Eq.  (4).  Calculations  of  the  trans- 
missivity as  a function  of  the  cirrus  cloud  thickness 
and  ice  content  were  carried  out  by  Liou  and  Stoffel 
(1976).  It  was  assumed  that  the  cirrus  cloud  was 
composed  of  long  circular  cylinders  randomly  oriented 
in  space  with  a base  height  of  8 km  and  an  isothermal 
temperature  of  — 36 °C.  Mean  wavenumbers  of  900, 
950,  11(K)  and  1150  cm-1  were  chosen  in  this  investiga- 
tion. Fig.  2 shows  the  ratios  of  the  transmissivities 
of  950,  1100  and  1150  cm-1  to  the  transmissivity  of 
900  cm-1.  It  is  evident  that  to  a good  approximation, 
these  curves  are  fairly  close  to  a linear  function  of 
the  thickness. 

On  the  basis  of  the  theoretical  radiative  transfer 
calculations,  we  postulate  that  the  transmissivities  of 
cirrus  clouds  for  wavenumbers  in  the  window  region 
can  be  scaled  as 

rj-  (a,As+6i)r*,  (10) 


Tablh  I.  Linear  coefficients  for  cloud  transmissivities. 


Hem'1) 

i 

a. 

bi 

900 

t 

0 

1 

950 

2 

-0.02 

0.99 

1100 

3 

-0.25 

1.23 

1150 

4 

-0.04 

1.01 

Ftc.  2.  Ratios  of  the  cloud  transmissivity  of  950,  1100  and 
1150  cm'1  to  that  of  900  cm"1.  The  cloud  transmissivities  are 
computed  from  the  radiative  transfer  program  including  scat- 
tering, absorption  and  emission  of  water  vapor  and  cloud 
particles. 


where  we  change  the  wavenumber  index  Ar  to  *. 
In  Eq.  (10),  tc  represents  the  transmissivity  at  a 
reference  wavenumber  (900  cm-1  in  Fig.  2),  and  a. 
and  bi  are  certain  constants  given  in  Table  1. 

In  view  of  the  linear  approximation  for  the  ratios 
of  the  transmissivities,  the  cloud  fraction  tj,  the  surface 
temperature  T„  the  cloud  thickness  As  and  the  cloud 
transmissivity  t‘  at  the  reference  wavenumber  may 
be  determined  from  four  radiance  measurements  in 
the  window  region  provided  that  the  atmospheric 
temperature  and  humidity  profiles  are  given,  for 
example,  from  climatological  data.  For  the  convenience 
of  the  following  discussions,  let 


/*<=  f fl,[7'(z)l/r,(zk,z)^ 

J 1-0 

r*=*r,(s»,0) 


(ID 


with  * denoting  values  to  be  evaluated  from  the  cli- 
matological data.  Thus  from  Eq.  (9)  we  have  the 
following  expressions  for  the  four  upwelling  radiances: 


7T“  [1  — n+>i(aiAs+fti)TeX5i(rt)r*+C]. 

*“1,2, 3, 4. 


Defining 


(?i(r.)-7r/[Bi(r.)r:+/:]. 


(12) 

(13) 
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Eq.  (12)  can  be  written  in  the  form 

1— 6.)rc,  .-1,2,3,  4,  (14) 

a nonlinear  equation  consisting  of  four  unknown  pa- 
rameters, 7,  t‘,  Az  anil  7',.  First,  we  tinil  that  n can 
be  eliminated  to  yield 

((>i  - 1 )[  (tf.A*+ft.)r‘  - 1 ] - ((>.  - 1 ) (r‘  - 1 ) = 0, 

» = 2,3, 4.  (15) 


Second,  r*  can  be  omitted  in  Eq.  (15)  to  obtain 

«?i-(W[0>i-l)(«.As+».  )-((>,-!)] 

-1=0, 

((>. -QdKQi  ~ 1 ) (asAz+M  - ((>, - 1 )] 

i = 3,4.  (16) 

Finally,  from  the  two  equations  in  (16)  we  eliminate 
Az  to  give 


[(>.(*.-*,•)+(>,(! (A 1 )]C(>i(«z -«4)-f-<?^4 -<  V/,] 

1=0.  (17) 

C(MN  — A;)+(L(1  l)][(h(<J.. — <J  ])  + ().•<! 


Eq.  (17)  represents  a nonlinear  equation  consisting 
of  an  unknown  parameter  T„  which  can  be  deter- 
mined from  a set  of  four  upwelling  radiances  /f‘  mea- 
sured in  a partially  cloudy  atmosphere.  On  deriving 
the  surface  temperature,  the  cloud  thickness,  the  cloud 
transmissivity  at  the  reference  wavenumber  and  the 
fraction  of  cloudiness  can  he  subsequently  evaluated. 
In  the  next  section,  we  shall  investigate  the  possible 
errors  of  the  resulting  parameters  caused  by  the  use 
of  the  climatological  temperature  and  humidity 
profiles  and  the  possible  random  errors  in  the  expected 
upwelling  radiances. 

3.  Computational  analyses 

In  order  to  test  the  idea  outlined  in  the  previous 
section,  computational  analyses  have  been  carried  out. 
Water  vapor  absorption  properties  are  first  discussed 
in  Section  3a.  In  3b,  a synthetic  partially  cloudy 
atmosphere  is  constructed.  Error  analyses  are  then 
presented  in  3c.  Some  notes  are  finally  given  in  3d 
for  the  estimation  of  the  vertical  ice  content  of  cir- 
rus clouds. 

a.  Water  taper  absorption  in  the  10  nm  window 

The  transmission  function  for  water  vapor  in  the 
10  jim  window  region  is  due  to  the  selective  absorption 
by  weak  lines  and  the  continuous  absorption,  and  can 
be  written 

r,  = r, (selective)  Xr, (continuum).  (18) 

The  selective  absorption  by  weak  lines  is  well  under- 
stood (see,  e.g.,  Goody,  1064).  However,  continuous 
absorption  in  the  10  jim  window  has  been  a subject 
of  considerable  speculation.  In  recent  years,  it  seems 
to  have  been  established  that  the  water  vapor  con- 
tinuum is  caused  by  the  broadening  of  foreign  gases 
and  the  e-type  (water  dimmer)  absorption  (Bignell, 
1070;  Burch,  1070).  Thus  the  absorption  coefficient  of 
the  continuum  may  be  expressed  as 

*(7>,e)  = *1(7')/>+*»(7>,  (19) 


where  p is  the  total  pressure  and  c the  partial  vapor 
pressure.  Based  on  BignelFs  measurements,  k:  is  much 
greater  than  k i,  although  the  temperature  dependence 
of  k ; has  not  been  investigated  completely. 

With  these  sources  of  absorption  taken  into  con- 
sideration the  transmission  function  for  water  vapor 
may  be  written  as 

r,(#)  = exp[— — (kip-\-ktf)u],  (20) 

where  u is  the  vertical  path  length  of  water  vapor 
and  the  weak  line  parameter 

y,=  (tou)*/8,  (21) 

with  o the  mean  line  intensity,  S the  mean  line  spacing 
and  a the  mean  line  width.  These  weak  line  parame- 
ters have  been  given  by  Roach  and  Goody  (1658) 
based  on  titling  the  statistical  band  model  from  ob- 
served data.  According  to  Bignell  (1970)  and  Burch 
(1970),  the  e-type  absorption  continuum  is  dominant 
in  the  10  jmi  window  region.  Furthermore,  compari- 
sons of  synthetic  IRIS  (Infrared  Interferometer  Spec- 
trometer) radiances  computed  from  Rigncll's  water 
vapor  continuum  data  with  observed  radiances  gave 
good  agreement  (Kunde  el  ill.,  1974).  In  view  of  this 
evidence  and  in  view  of  the  fact  that  the  absorption 
coefficient  is  continuous,  we  have  computed  water 
vapor  transmissivities  employing  the  mean  wave- 
numbers  of  900,  950,  1100  and  1150  cm  1 in  con- 
junction with  radiative  transfer  calculations  reported 
in  this  paper. 

b.  Synthetic  partially  cloudy  atmosphere 

To  test  the  accuracy  of  the  procedures  described 
in  Section  2,  a synthetic  partially  cloudy  atmosphere 
lias  been  constructed  (Fig.  3).  It  is  assumed  that 
within  the  field  of  view  of  the  satellite  radiometer, 
there  is  70^)>  of  cirrus  cloudiness  whose  thickness  is 
2 km,  and  the  base  of  the  cirrus  cloud  is  located  at 
8 km  in  the  atmosphere.  Based  on  radiative  transfer 
calculations,  we  obtain  a transmissivity  of  0.675  for 
a 2 km  cirrus  cloud  at  a wavenumber  of  900  cm  '. 
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Moreover,  a surface  temperature  of  2V4  K is  employed. 
With  res|H-it  to  the  atmospheric  ten\|>erature  and 
humidity  proldes,  we  use  the  tahulated  climatological 
data  l>t  Mcl'lalchev  el  a/.  (1V71)  for  the  mid-latitude 
summer  atmosphere. 

With  the  transmissivity  at  ‘Ml  cm  1 given,  the  rest 
of  the  transmissivities  can  he  evaluated  according  to 
Kq  (It"  lleme,  we  can  theoretically  calculate  the 
expected  upwelling  radiances  in  the  partially  cloudy 
atmosphere  for  four  wavenumbers  !>\  means  of  Kqs. 
(I'  (.!'  Ome  the  upwelling  radiances  have  been  ob- 
tained, we  proceed  to  Kqs.  (14)  (17'  to  retrieve  the 
surface  teni|N-rature,  the  cloud  thickness,  the  trans- 
missivity at  VtKi  cm  1 and  the  fraction  of  cirrus 
cloudiness. 

r.  h'rror  analyses 

F.q  (171  represents  a complicated  nonlinear  equa- 
tion which  contains  the  unknown  surface  temperature. 
We  may  write 

b\r,)-0.  (22) 

Since  the  function  is  a high  order  nonlinear  equation, 
we  would  e\|iect  that  it  consists  of  a number  of  roots 
for  f,  in  the  (tl, * ) domain.  However,  there  will  he 
only  one  root  that  is  our  desired  solution,  big.  4 
illustrates  the  behavior  of  F.q.  (221  for  the  conditions 
given  in  Seition  .lb  with  the  realistic  range  of 


from  240  to  ,140  K.  In  addition  to  the  desired  root 
of  2V4  K,  we  find  that  2b. 1 K also  satisties  Kq.  (22). 
Hence,  care  should  be  taken  in  the  error  analysis  to 
disregaril  the  undesirable  roots.  Fortunately,  from  a 
number  of  error  analyses  exercises,  we  found  that  the 
correct  value  always  lies  at  the  xero  intersection  of 
the  right-hand  curve  in  Fig.  4. 

The  first  error  analyses  were  done  for  the  upwelling 
radiance.  Random  numbers  were  arbitrarily  selected 
and  added  to  the  upwelling  radiances  calculated 
exactly  from  the  synthetic  atmosphere.  The  retrieval 
procedures  using  the  upwelling  radiances  with  random 
errors  were  carried  out  for  the  surface  temperature, 
and  subsequently  for  the  cloud  thickness,  the  reference 
cloud  transmissivity  and  the  cloud  fraction.  We  define 
the  following  statistical  terms: 

1 V 

A/,„—  £ A/.//..  (2.1) 

X 

SI -[£  (24) 

i-i 

Table  2 lists  the  resulting  parameters  in  terms  of 
these  two  variables.  The  second  line  of  this  table 
shows  that  1.2'-,  random  error  produces  almost  negli 
gible  errors  in  the  four  resulting  parameters.  A random 
error  of  11. 2*' J.  introduces  deviations  of  about  .l.b  k, 
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A/„. 

T4°  k) 

A/(Wm) 

t 

rINNI 

a 

tl 

0 

2**4 

i 

0.075 

0.7 

0.012 

0 002 

204  .140 

1 024 

0 084 

0 702 

0 0So 

0.010 

70S  745 

1 0S5 

0 710 

0.71(1 

0 112 

0.020 

207  S40 

1 385 

0 760 

0.710 

(1002 

0 003 

207  480 

2 187 

0 044 

0 60(1 

t>.fi  1 km,  0.08  itntl  0.02  for  tin'  surface  temperature, 
the  cloud  thickness,  the  reference  cloud  transmissivity 
and  the  fraction  of  cloudiness,  respectively.  Other 
experiments  employing  different  sets  of  random  num- 
bers were  also  carried  out.  It  was  found  that  a large 
value  of  A/rm,  would  introduce  large  errors  in  the 
four  resulting  parameters  The  fifth  line  of  this  table 
illustrates  an  example  as  to  how  the  root  mean  square 
affects  the  present  analysis.  We  see  that  a small  random 
error  of  0.2*',  having  a larger  mis  (l)J%)  produces 
about  the  same  errors  in  the  four  resulting  parameters 
as  the  case  when  a random  error  of  11.2%  with  a 
small  rms  of  2%  is  introduced. 

It  is  known  that  the  temperature  and  humidity 
profiles  derived  from  radiosonde  are  normally  not 
available  in  ocean  areas,  particularly  in  the  Southern 
Hemisphere.  The  best  guesses  of  the  temperature  anil 
humidity  profiles  in  cloudy  conditions,  perhaps,  may 
be  obtained  from  climatology.  Deviations  of  the 
climatological  values  from  the  true  profiles  are  likely 
to  occur.  Thus,  it  is  of  importance  to  examine  how 
the  errors  in  the  temperature  and  humidity  profiles 
affect  the  resulting  recovered  parameters  develojred 
previously.  (Note  that  if  true  temperature  and  hu- 


midity profiles  are  available,  such  error  analyses  are 
obviously  not  needed.) 

In  order  to  investigate  the  effect  of  utilizing  dima 
tological  temperature  and  humidity  profiles,  three 
experiments  were  carried  out:  (i)  random  errors  from 
(I  to  5 k and  from  t>  to  It)  k were  added  to  the  tem- 
perature profile,  (ii > random  errors  from  tl  to  5% 
and  from  (I  to  UK’,"  were  added  to  the  relative  hu- 
midity profile,  and  (iii)  independent  random  errors 
similar  to  (i)  and  (ii)  for  both  temperature  and  hu- 
midity profiles  were  used.  Results  are  shown  in  f ig  5. 
terrors  in  the  temperature  profile  apparently  introduce 
large  errors  in  the  resulting  four  parameters.  Kor 
random  errors  from  0 to  5 k in  the  temperature 
profile,  we  note  that  deviations  of  about  V k,  1 km, 
0.07  and  O.tMKi  are  produced  for  the  surface  tempera 
lure,  the  cloud  thickness,  the  reference  Iransmissivitx 
and  the  cloud  fraction,  respectively.  On  the  other 
hand,  random  errors  from  tt  to  s',  in  the  Immiditv 
profiles  introduce  s k.  O s km,  (Mil  and  tl.dtl.t  for  the 
resulting  four  parameters,  respectively  . Note  that  the 
+ and  — signs  in  these  diagrams  indicate  the  positive 
and  negative  errors,  respectively . When  independent 
errors  are  added  to  both  the  temperature  and  hu- 
midity protiles,  we  find  that  errors  in  the  four  resulting 
parameters  are  greatly  reduced.  The  reason  for  this 
may  be  caused  by  the  compensating  effect  of  opposite 
errors  for  the  temperature  and  humidity  profiles  as 
evident  from  fig.  5.  It  should  lie  noted  that  simul- 
taneous deviations  of  the  true  tenqierature  and  hu- 
midity profiles  from  the  climatological  means  are  to 
be  expected  in  the  atmosphere.  Hence  we  would 
anticipate  that  the  use  of  climatological  data  in  the 


Fie,.  4.  Behavior  of  the  nonlinear  equation  (17)  as  a function  of 
the  surface  temperature. 
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present  titiuv.il  analysis  iutriHltu  vs  iusigiiihiaiil  errors 
m i limit  jiaraiiietcrs. 

In  tliv  development  t>(  tin-  rftricv.il  technique,  the 
assumption  h.is  I uiii  m.utf  tli.it  the  cloud  is  hi^li  m 
the  atmosphere  vi  tli.it  radiative  elicits  of  water  vapor 
.iImivc  the  ilmiil  i .in  In-  neglected  Heine  the  method 
should  In  .i|i|ihisl  ,is  sin  h and  i uses  w hen  Inn  cumulus 
.uni  nr  str.itus  .ire  present  .ue  to  lie  disregarded. 
Suin'  visthle  .uni  IK  ilmnl  picture*  .ire  routinely 
.iv.ul.ililf  Inmi  nperAtiuii.il  satellites,  iltsi  i iunn.it inn 
In-tween  Inw  .uni  high  ilnuils  appeals  possible  1‘liere 
are  uses  when  high  ami  Inw  ilnuils  occur  siniulta 
nemislx  in  the  atmosphere  In  these  situations,  low 
ilnuils.  whnli  are  generally  optnalli  thiih,  max  he 
iniisiiletiil  as  suilans  ami  to  a good  approximation 
the  theorelii.il  approach  rc|>orled  here  ma\  also  be 
appln  ahh  Note  that  we  have  eniploM-il  'MHI,  ‘)s|l, 
Him  ami  n>t'  i m 1 in  the  theoretif.il  analysis  to 
ilcmniistiati  that  ililormalion  associated  with  iiims 
mai  he  dimed  limn  satellite  i.uliame  nltserv.it iolls, 
t M loutsf,  otliei  spfitral  hamls  within  the  wimlow 
I if  ion  oufht  to  In-  explored  to  111x1111-  I lie  optimum 
hamls,  wliii  h could  In  pioposeil  loi  iuture  satellite 
experiments  cniiierning  the  ilelctlinl)  ol  urrus  cloud 
composition*  ami  structure  in  the  atmosphere. 

1 /.  .Vt>/e  on  ili<  ilttffmiHiitum  oj  urm.i  it e uwtent 

In  the  preceding  anahses,  we  have  shown  that  the 
surface  temperature,  the  cloud  thickness,  the  reference 
ilmnl  transmissivity  ami  the  fraction  of  cloudiness 
mac  he  inferred  from  four  radiance  measurements  in 
the  window  region.  Once  a reliable  reference  cloud 
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transmissivity  has  been  denied,  it  seems  possible  to 
estimate  the  vertical  in-  iiiiiteiil  following  the  proce- 
dures described  below 

Itascd  on  the  theou-tical  radiative  transfer  1 ah  illa- 
tions, we  max  construct  a diagram  fl  irt  <>)  consisting 
ot  the  transmissivity  at  a given  wavenumber  and  the 
optical  depth  at  that  ftvcii  wavetiumhei  We  note 
that  the  transmissivity  is  lint hiii^*  but  an  attenuation 
parameter  flor  a urrus  ilouit  whose  optical  depth 
is  less  than  about  I'  ' it  appears  that  an  exponential 
attenuation  max  he  appln-il  ' |(\  means  ot  this  curve 
1 oust rin  ted  prior  to  tin-  retrieval  experiment,  a trails 
nussivitx  derived  from  the  previous  analysis  should 
correspond  to  an  optical  depth 

We  now  express  the  optical  depth  as 

r(H-A!d«i(H.  (25) 

where  d.,,(i'  denotes  the  volume  cxlimtioti  cross 
sci  lion  tor  a given  mean  wavenumber  1 Now  with 
the  information  ot  the  cloud  thickness  wliii h has  been 
determined  simultaneously  with  the  1 loud  Ir.insmis 
six  1 1 > . the  important  volume  extinction  uoss  section 
may  be  estimated. 

Since  cirrus  clouds  are  composed  ice  parities  on 
the  order  of  ItXI  IlHX'  >»m  (see,  eg,  llevmslield  and 
knollcnberg,  1072),  an  extinction  parameter  (usually 
denoted  as  ( *,.» 1 1 of  about  2 may  be  adopted  Con- 
sequently, the  extinction  cross  section  is  mainly  a 
function  of  the  size  and  concentration  of  ice  particles 
and  may  be  written  as 

d«i~()«i-t.V,  (2ft) 


Km  V t-'rrur  analyses  fur  the  use  ut  the  dimatwlwgical  temperature  anil  humidity  |irutilea. 
The  almisaa  anil  animate  represent  raiuluni  etixirs  amt  deviations  uf  the  lour  parameters 
Irom  the  true  values,  respectively. 
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Fir*,  ft.  Transmissivity  as  a function  of  optical  depth  Utsed  on 
radiative  transfer  calculations. 
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where  X anil  .{  denote  the  concentration  anil  the 
mean  cross  section  area  of  ice  particles,  respectively. 
For  ice  crystals  randomly  oriented  in  a horizontal 
plane,  we  should  have  ,i~/i7,  with  7 and  <7  represent- 
ing the  mean  length  in  the  major  and  minor  axes, 
respectively. 

Finally,  the  vertical  ice  content  can  be  written  in 
the  form 

K' *pif\YAs,  (27) 

where  pi  represents  the  density  of  the  ice  and  f‘  the 
mean  volume  of  ice  particles. 

If  we  let  $ F .1,  the  above  equation  for  the  vertical 
ice  content  is  given  by 

0™, 

ll'»  pi fAs.  (28) 

(>„, 

For  long  circular  ice  cylinders  randomly  oriented  in 
a horizontal  plane,  we  have  £ = «7  4.  Here,  pi,  (h,i 
and  £ are  known  quantities.  Thus,  a determination 
of  and  A:  gives  an  estimation  on  the  vertical 
ice  content. 

4.  Conclusion 

We  have  presented  a retrieval  method  for  the  de- 
termination of  surface  temperature,  cirrus  cloud  thick- 
ness, cirrus  cloud  transmissivity  at  a reference  wave 
number  and  the  fraction  of  cirrus  cloudiness  within 
the  held  of  view  from  four  cxpeitcd  upwelling  radiance 
observations  in  the  |t>  pm  window  region.  Theoretical 
analyses  assume  that  water  vapor  effects  above  the 


cirrus  cloud  are  negligible  and  that  the  ratios  of 
transmissivities  are  linear  functions  of  the  cloud 
thickness.  Justifications  of  these  two  assumptions  are 
given  on  the  basis  of  radiative  transfer  calculations. 

Computational  analyses  employing  the  climatological 
means  show  that  random  errors  in  both  the  tem- 
perature and  humidity  profiles  produce  insignificant 
errors  in  the  four  resulting  parameters.  Random  errors 
in  the  expected  upwelling  radiances  also  introduce 
small  errors  in  the  resulting  parameters  provided  that 
I lie  rms  random  errors  are  relatively  small  as  com- 
pared with  the  mean  value.  Hence,  from  a theoretical 
point  of  view,  the  retrieval  procedures  demonstrated 
appear  feasible  for  the  evaluation  of  cirrus  cloud 
thickness  and  transmissivity  at  a reference  wave- 
number  in  the  10  pm  window.  Once  the  thickness 
and  transmissivity  at  a given  wavenumber  of  a cirrus 
cloud  have  been  derived,  we  note  that  the  vertical 
ice  content  of  that  cirrus  cloud  may  he  estimated. 

The  retrieval  technique  described  in  the  preceding 
sections,  in  principle,  should  be  applicable  to  cirrus 
cloudy  atmospheres  for  the  quantitative  estimation 
of  cirrus  thickness  and  ice  content  from  passive  satel- 
lite observations.  Although  a ground-based  active 
remote  sensing  method  employing  a laser  source  has 
been  recently  developed  for  the  detection  of  cirrus, 
it  is  of  local  nature  and  is  to  be  verified  for  the  suc- 
cessful determination  of  thickness  and  ice  content. 
However,  cirrus  clouds  cover  a large  portion  of  the 
planet.  It  is  therefore  vitally  important  to  explore 
passive  remote  sensing  techniques  from  which  cirrus 
cloud  composition  and  structure  may  be  derived  front 
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satellites  on  a routine  busts  over  the  global  settle. 
There  has  been  hardly  any  investigation  on  the 
recovery  of  cloud  properties  front  passive  satellite 
sensing  owing  to  the  complexity  of  the  cloud  inter- 
action with  the  radiation  held  of  the  atmosphere. 
This  is  particular!)  apparent  for  the  high,  semi-trans- 
parent lirrus  clouds.  This  palter  represents  an  approach 
based  on  which  the  thickness  ami  ice  content  of  cirrus 
may  be  evaluated  by  utilizing  four  pieces  of  radiance 
information  in  the  10*<m  window  region. 

Any  satellite  sensing  technique  for  the  recovery  of 
atmospheric  composition  and  structure  is  subject  to 
verification  by  utilizing  real  data.  The  data  source 
that  may  lie  of  use  in  connection  with  the  present 
theoretical  analysis  would  be  the  IRIS  experiment 
on  board  Nimbus  4 as  reported  by  kunde  ft  ill.  (l*)74'. 
IRIS  measured  the  thermal  emission  of  the  earth's 
atmosphere  and  surface  from  400  1600  cur1  with  an 
a|>odized  spectral  resolution  of  2.8  cm-1.  It  seems 
possible  to  select  a set  of  four  spectral  intervals  in 
the  750  1250  cm  1 region  to  earn  out  the  analyses 
described  in  this  paper  provided  that  the  synoptic 
and  cloud  information  of  a cirrus  cloudy  atmosphere 
is  available.  I nforlunatelv,  satellite  cloud  experiments 
require  in  situ  information  of  the  cloud  thickness  and 
composition,  which  are  normally  not  available  under 
the  satellite  pass.  Perhaps  the  major  thrust  of  the 
present  paper  is  to  illustrate  that  certain  cirrus  cloud 
parameters  could  be  evaluated  front  a set  of  radiance 
observations  in  the  10  pm  window  region.  Such  a 
theoretical  demonstration  seems  to  warrant  a pos- 
sibility of  pre-satellite  ex|>eriments  for  the  quanti- 
tative determination  of  cirrus  cloud  composition  and 
thickness. 
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